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ON  THE  GORTLER  INSTABILITY  IN  HYPERSONIC  FLOWS: 
SUTHERLAND  LAW  FLUIDS  AND  REAL  GAS  EFFECTS  1 


Yibin  B.  Fu,  Philip  Hall,  and  Nicholas  D.  Blackaby 
Department  of  Mathematics 
Oxford  Road 
University  of  Manchester 
Manchester  M13  9PL 
U.K. 

Abstract 

The  Gortler  vortex  instability  mechanism  in  a  hypersonic  boundary  layer  on  a  curved  wall 
is  investigated  in  this  paper.  Our  aim  is  to  clarify  the  precise  roles  of  the  effects  of  boundary 
layer  growth,  wall  cooling  and  gas  dissociation  in  the  determination  of  stability  properties. 
We  first  assume  that  the  fluid  is  an  ideal  gas  with  viscosity  given  by  Sutherland’s  law.  It 
is  shown  that  when  the  free  stream  Mach  number  M  is  large,  the  boundary  layer  divides 
into  two  sublayers:  a  wall  layer  of  0(M 3^2)  thickness  over  which  the  basic  state  temperature 
is  0(M2)  and  a  temperature  adjustment  layer  of  0(1)  thickness  over  which  the  basic  state 
temperature  decreases  monotonically  to  its  free  stream  value.  Gortler  vortices  which  have 
wavelength  comparable  with  the  boundary  layer  thickness  (i.e.  have  local  wavenumber  of 
order  M~3/2)  are  referred  to  as  wall  modes.  We  show  that  their  downstream  evolution  is 
governed  by  a  set  of  parabolic  partial  differential  equations  and  that  they  have  the  usual 
features  of  Gortler  vortices  in  incompressible  boundary  layers.  As  the  local  wavenumber 
increases,  the  neutral  Gortler  number  decreases  and  the  centre  of  vortex  activity  moves 
towards  the  temperature  adjustment  layer.  Gortler  vortices  with  wavenumber  of  order  one 
or  larger  must  necessarily  be  trapped  in  the  temperature  adjustment  layer  and  it  is  this 
mode  which  is  the  most  dangerous.  For  this  mode,  we  find  that  the  leading  order  term  in 
the  Gortler  number  expansion  is  independent  of  the  wavenumber  and  is  due  to  the  curvature 
of  the  basic  state.  This  term  is  also  the  asymptotic  limit  of  the  neutral  Gortler  numbers 
of  the  wall  mode.  To  determine  the  higher  order  correction  terms  in  the  Gortler  number 
expansion,  we  have  to  distinguish  between  two  wall  curvature  cases.  When  the  wall  curvature 
is  proportional  to  (2x)~3^2  where  x  is  the  streamwise  variable,  the  Mach  number  M  can  be 
scaled  out  of  the  problem  and  we  show  that  in  the  0(1)  wavenumber  regime,  Gortler  vortices 

'This  research  was  supported  in  part  by  the  National  Aeronautics  and  Space  Administration  under  NASA 
Contract  No.  NASI- 18605  while  the  author  was  in  residence  at  the  Institute  for  Computer  Applications  in 
Science  and  Engineering  (ICASE),  NASA  Langley  Research  Center,  Hampton,  VA  23665. Additional  support 
was  provided  by  USAF  under  Grant  AFOSR89-0042  and  SERC. 
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are  again  governed  by  a  set  of  parabolic  partial  differential  equations  and  therefore  the  higher 
order  correction  terms  in  the  Gortler  number  expansion  are  not  uniquely  determined  and 
are  strongly  dependent  on  nonparallel  effects.  In  the  large  wavenumber  limit,  however, 
nonparallel  effects  become  of  second  order;  Gortler  vortices  evolve  downstream  in  a  quasi¬ 
parallel  manner  and  the  Gortler  number  expansion  has  its  first  three  terms  independent  of 
nonparallel  effects.  In  the  more  general  case  when  the  wall  curvature  is  not  proportional  to 
(2a:)-3/2,  the  effect  of  the  curvature  of  the  basic  state  persists  in  the  downstream  development 
of  Gortler  vortices;  non-parallel  effects  are  important  over  a  larger  range  of  wavenumbers  and 
they  become  of  second  order  only  when  the  wavenumber  is  of  order  higher  than  0(M 1/4). 
In  the  latter  case  the  Gortler  number  expansion  has  the  first  two  terms  independent  of 
nonparallel  effects;  the  first  term  being  due  to  the  curvature  of  the  basic  state  and  the 
second  term  due  to  viscous  effects.  The  second  term  becomes  comparable  with  the  first  term 
when  the  wavenumber  reaches  the  order  M3,\  in  which  case  another  correction  term  can 
also  be  found  independently  of  nonparallel  effects.  Next  we  investigate  real  gas  effects  by 
assuming  that  the  fluid  is  an  ideal  dissociating  gas.  We  find  that  both  gas  dissociation  and 
wall  cooling  are  destabilizing  for  the  mode  trapped  in  the  temperature  adjustment  layer, 
but  for  the  wall  mode  trapped  near  the  wall  the  effect  of  gas  dissociation  can  be  either 
destabilizing  or  stabilizing. 
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1  Introduction 


This  paper  is  an  extension  of  our  previous  paper  Hall  and  Fu  (1989)  on  the  linear  development 
of  Gortler  vortices  and  the  reader  is  referred  to  that  paper  for  a  more  detailed  review  of  the 
relevant  literature.  In  that  paper,  we  assumed  that  the  fluid  was  an  ideal  gas  with  the  viscosity 
given  by  Chapman’s  law.  It  was  found  that  when  a  hypersonic  boundary  layer  first  loses 
stability  to  Gortler  vortices,  the  vortices  are  necessarily  trapped  in  the  logarithmically  thin 
temperature  adjustment  layer  over  which  the  temperature  of  the  basic  flow  changes  rapidly 
to  its  free  stream  value.  In  other  words,  the  mode  trapped  in  the  temperature  adjustment 
layer  has  a  smaller  Gortler  number  than  any  other  mode.  As  a  consequence  of  this  region  of 
vortex  activity  being  thin  (which  leads  us  to  consider  Gortler  vortices  of  small  wavelength), 
the  perturbation  equations  governing  the  downstream  development  of  the  vortices  reduce  to 
ordinary  differential  equations  within  the  order  of  approximation  considered  if  the  appropriate 
“fast”  streamwise  dependence  of  the  instability  is  built  into  the  disturbance  flow  structure. 
Thus  the  non-uniqueness  of  the  neutral  stability  curve  associated  with  incompressible  Gortler 
vortices  disappears  at  high  Mach  numbers  and  a  unique  neutral  curve  with  distinct  left  and 
right  branches  is  obtained. 

However,  a  real  fluid  has  its  viscosity  given  by  the  more  complicated  Sutherland’s  law. 
Although  in  most  of  the  previous  investigations  on  compressible  boundary  layers  Chapman’s 
viscosity  law  has  been  adopted  as  an  approximation  to  Sutherland’s  law,  such  an  approximation 
is  poor  for  hypersonic  flows  in  which  the  fluid  temperature  varies  significantly  across  the 
boundary  layer.  Thus  it  is  of  interest  to  investigate  how  our  previous  results  are  modified 
if  the  more  realistic  Sutherland’s  law  is  adopted.  This  is  one  of  the  problems  which  we  are 
addressing  in  the  present  paper. 

The  other  problems  which  we  consider  are  the  effects  of  gas  dissociation  and  wall  cooling 
on  the  flow  stability.  For  a  hypersonic  boundary  layer,  the  temperature  near  the  wall  is 
typically  of  order  0(M2)  where  M  is  the  free  stream  Mach  number,  and  gas  dissociation  must 
necessarily  take  place.  Also,  in  practical  situations,  walls  can  not  possibly  withstand  such  high 
temperatures  and  they  must  be  cooled.  Thus  it  is  also  of  special  interest  to  clarify  the  precise 
roles  of  these  two  mechanisms  in  the  stability  properties  of  hypersonic  boundary  layers. 

The  major  difference  between  Gortler  vortices  in  incompressible  and  hypersonic  flows  is  that 
the  presence  of  the  temperature  adjustment  layer  at  the  edge  of  a  hypersonic  boundary  layer 
where  the  basic  temperature  field  decreases  rapidly  to  its  free-stream  value  enables  hypersonic 
Gortler  vortices  to  be  concentrated  well  away  from  the  wall.  In  the  incompressible  case  we 
know  from  the  work  of  Hall  (1982,1983),  and  Denier,  Hall  and  Seddougui  (1990)  that  at  order 
one  Gortler  numbers  unstable  Gortler  vortices  are  not  localized  within  the  basic  boundary 
layer.  At  higher  Gortler  number  the  most  dangerous  Gortler  vortices  have  a  wavelength  small 
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compared  to  the  boundary  layer  thickness  and  are  trapped  near  the  wall.  At  order  one  Mach 
numbers  this  situation  does  not  change  significantly  and  the  nonparallel  problem  has  been 
discussed  by  Wadey  (1990)  and  Spall  and  Malik  (1989).  In  the  latter  two  investigations  the 
nonparallel  equations  were  solved  numerically  following  the  approach  of  Hall  (1983),  the  main 
result  obtained  was  that  the  neutral  position  or  growth  rate  of  a  Gortler  vortex  is  a  function 
of  it’s  upstream  history.  However  the  numerical  calculations  of  Wadey  (1990)  suggest  that  as 
the  Mach  number  increases  the  position  where  an  unstable  Gortler  vortex  locates  itself  moves 
towards  the  edge  of  the  boundary  layer.  That  result  is  consistent  with  what  we  shall  find  in 
this  paper. 

The  present  paper  is  limited  only  to  the  linear  regime  of  vortex  growth  ;  nonlinear  aspects  of 
incompressible  or  low  Mach  number  Gortler  vortex  growth  have  been  discussed  by  Hall  (1988), 
and  Wadey  (1990).  For  a  detailed  account  of  the  nonlinear  regime  the  reader  is  referred  to 
the  review  paper  by  Hall  (1990).  Before  going  on  to  discuss  the  work  presented  here  we 
also  note  that  hypersonic  boundary  layers  are  also  susceptible  to  instabilities  not  induced  by 
streamline  curvature.  Thus,  for  example,  Cowley  and  Hail  (1990),  Blackaby,  Cowley  and  Hall 
(1990),  Smith  and  Brown  (1990)  have  discussed  the  role  of  Rayleigh  or  Tollmien-Schlichting 
wave  instabilities  in  hypersonic  boundary  layers.  Thus  any  nonlinear  investigation  of  Gortler 
vortices  at  hypersonic  speeds  must  allow  for  the  possible  interaction  of  the  vortices  with  other 
finite  amplitude  instability  mechanisms. 

This  paper  is  organized  as  follows.  In  §2  and  in  the  first  part  of  §3  we  discuss  the  basic 
state.  Here  we  discuss  the  significant  changes  in  the  basic  state  which  occur  when  Sutherland’s 
law  is  used  instead  of  the  Chapman  law.  In  particular,  the  logarithmically  thin  temperature 
adjustment  layer  found  for  Chapman  fluids  is  now  replaced  by  a  more  complex  adjustment  layer 
of  0(1)  thickness.  We  shall  see  later  that  this  difference  will  strongly  affect  the  downstream 
evolution  properties  of  Gortler  vortices. 

In  §3  we  formulate  the  linear  stability  problem  for  a  hypersonic  boundary  layer  and  the 
linear  perturbation  equations  are  obtained  in  the  usual  manner  by  superimposing  a  Gortler 
vortex  structure  on  the  basic  state  and  linearizing  the  Navier-Stokes  equations.  When  Chap¬ 
man’s  law  is  used,  these  perturbation  equations  can  be  reduced  to  a  set  of  ordinary  differential 
equations  after  the  wavelength  of  the  vortices  is  scaled  by  the  thickness  of  the  logarithmically 
thin  adjustment  layer.  However,  here  we  show  that  because  the  temperature  adjustment  layer 
is  of  0(1)  thickness  when  Sutherland’s  law  is  used,  non-parallel  effects  are  more  pronounced 
and  their  effects  are  different  for  different  wall  curvatures.  To  be  more  specific,  when  the 
wall  curvature  is  proportional  to  ( 2x)~ 3/<2,  the  curvature  of  the  basic  state  only  gives  rise  to 
an  0(Af3/2)  wavenumber  independent  term  in  the  Gortler  number  expansion  and  its  effect  is 
not  present  in  the  downstream  development  of  Gortler  vortices  in  the  neighbourhood  of  the 
neutral  position;  whilst  in  the  more  general  case  when  the  wall  curvature  is  not  proportional 
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to  (2x)“3/2,  the  curvature  of  the  basic  state  not  only  gives  rise  to  an  0(M3^2)  wavenumber  in¬ 
dependent  term  in  the  Gortler  number  expansion,  but  it  also  affects  the  downstream  evolution 
of  Gortler  vortices  in  the  neighbourhood  of  the  neutral  position  and  thus  affects  the  determi¬ 
nation  of  other  higher  order  correction  terms  to  the  Gortler  number  expansion.  Sections  4  and 
5  are  respectively  devoted  to  the  discussion  of  these  two  cases. 

In  section  4,  we  show  that  in  the  special  curvature  case  /c(x)  =  (2x)~3/2,  the  Mach  number 
can  be  scaled  out  of  the  problem;  in  the  0(1)  wavenumber  regime,  the  perturbation  equa¬ 
tions  are  partial  differential  equations  and  they  have  to  be  solved  numerically  by  a  marching 
procedure.  We  present  our  numerical  results  which  show  that  neutral  curves  depend  crucially 
on  what  initial  perturbations  we  impose  and  where  we  impose  them.  In  the  large  wavenum¬ 
ber  limit,  nonparallel  effects  are  negligible  and  a  simple  asymptotic  expression  is  obtained  for 
the  Gortler  number  in  terms  of  the  wavenumber.  In  section  5,  we  show  that  because  of  the 
persistent  effect  of  the  curvature  of  the  basic  state,  non-parallel  effects  are  important  over  a 
larger  range  of  wavenumbers  and  they  become  negligible  only  when  the  wavenumber  is  of  order 
larger  than  0(M 1/4).  Thus  for  wavenumbers  of  order  0(Ma)  with  a  <  1/4,  the  perturbation 
equations  which  govern  the  downstream  evolution  of  Gortler  vortices  are  partial  differential 
equations  and  the  situation  is  similar  to  the  0(1)  wavenumber  case  discussed  in  section  4. 
When  a  >  1/4,  nonparallel  effects  are  not  so  pronounced  and  the  Gortler  number  expansion 
has  the  first  two  terms  independent  of  nonparallel  effects;  the  first  term  due  to  the  curvature 
of  the  basic  state  and  the  second  term  due  to  viscous  effects.  The  second  term  becomes  com¬ 
parable  with  the  first  term  when  the  wavenumber  reaches  the  order  0(M 3/8),  in  which  case 
another  correction  term  can  also  be  found  independently  of  nonparallel  effects. 

To  complete  our  stability  analysis,  we  devote  section  6  to  the  wall  mode  which  has  wave¬ 
length  comparable  with  the  boundary  layer  thickness.  This  mode  is  nonparallel  and  neutral 
curves  have  to  be  obtained  by  solving  a  set  of  partial  differential  equations.  Our  numerical 
results  show  that  neutral  curves,  although  non-unique,  all  decrease  monotonically  with  the 
wavenumber  and  tend  to  a  constant  value  in  the  large  wavenumber  limit,  thus  matching  in  the 
large  wavenumber  limit  with  the  mode  trapped  in  the  temperature  adjustment  layer. 

In  section  7,  we  investigate  real  gas  effects  and  wall  cooling  effects.  We  assume  that  the 
fluid  is  an  ideal  dissociating  gas.  After  dissociation  has  taken  place,  the  fluid  becomes  a  gas 
mixture.  We  first  determine  the  constitutive  properties  of  the  gas  mixture  and  then  show  how 
our  previous  results  for  ideal  gases  are  modified  when  gas  dissociation  is  taken  into  account. 
We  show  for  the  mode  trapped  in  the  temperature  adjustment  layer  that  the  leading  order 
Gortler  number  is  decreased  by  both  gas  dissociation  and  wall  cooling  and  thus  we  conclude 
that  both  these  mechanisms  are  destabilizing.  For  the  wall  mode,  neutral  curves  are  not 
unique  and  so  we  cannot  draw  any  general  conclusion.  For  the  case  we  consider,  the  neutral 
curves  corresponding  to  the  two  models  intersect,  so  the  effect  of  gas  dissociation  can  be  either 
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destabilizing  or  destabilizing.  Finally  in  the  last  section  we  give  some  further  discussion. 


2  Basic  state 


Consider  a  hypersonic  boundary  layer  over  a  rigid  wall  of  variable  curvature  (l/A )k(x*/L), 
where  L  is  a  typical  streamwise  length  scale  and  A  is  a  lengthscale  characterizing  the  radius  of 
curvature  of  the  wall.  We  choose  a  curvilinear  coordinate  system  (x*,y*>z*)  with  x*  measuring 
distance  along  the  wall,  y*  perpendicular  to  the  wall  and  z*  in  the  spanwise  direction.  The 
corresponding  velocity  components  are  denoted  by  (u*,v*,w*)  and  density,  temperature  and 
viscosity  by  p*,T*  and  p*  respectively.  The  free  stream  values  of  these  quantities  will  be 
signified  by  a  subscript  oo.  We  define  a  curvature  parameter  6  by 


6  = 


L 

A' 


(2.1) 


and  consider  the  limit  6  -*  0  with  the  Reynolds  number  R  defined  by 

R  _  <o  LPlo 
Vlo 

taken  to  be  large  so  that  the  Gortler  number 


(2.2) 


G  =  2  R1/26 


(2.3) 


is  0(1).  In  the  following  analysis,  coordinates  (x*,y*,z*)  are  scaled  on  (£*,  R~l/2L,  R~1/2L), 
the  velocity  (u*,v*,w*)  is  scaled  on  (u^,  R~1^2ul0,  R~1/2ula)  and  other  quantities  such  as 
p*,T*,  and  p*  are  scaled  on  their  free  stream  values  with  the  only  exception  that  the  pressure  p* 
is  scaled  on  p^ouoo  and  the  coefficient  of  heat  conduction  k*  is  scaled  on  p*^.  All  dimensionless 
quantities  will  be  denoted  by  the  same  letters  without  a  superscript  *.  Then  the  Navier-Stokes 


equations  are  given  by 


dp  d  .  . 

at  +  =  °- 

Du  dp  d  .  du.  d  .  du. 
p~Dt  ~  ~~di  +  ~dy'P'di)  +  dx^dz^ 

"m  +  2Gku  =  -R%  +  ^  |(A  -  3"^)  +  top 

d  .  dv.  d  .  dv. 

+  dy^dy^  +  dz^dz^ 


dvp, 
dy  ' 


Dw 

~Dt 


dp  d  J.  2  .dvp)  d  .  dvp  d  .  dw  d  .  dw 
Redz  +  dz\(  ^dxpj+dxp^dz^dy^dy^dz^dz^ 


pcp 


DT  =  *  -  DM’[(g)’  +  (£}’]  +  (7  - 


Dt 


(2.4) 

(2.5) 


(2.6) 

(2.7) 


4 


I  8  ,.ars  Id  „8T. 

+  a  dy(k  dy^  +  a  dz^  dz^ 

7  M2p  —  (1  +  a)pT.  (2.9) 

Here  we  have  used  a  mixed  notation  in  which  (vi,U2,V3)  is  identified  with  ( u,v,w )  and 
(xj,  12,2:3)  with  ( x,y,z ).  Repeated  suffices  signify  summation  from  1  to  3.  The  functions 
X,k,cp  and  h  denote  in  turn  the  bulk  viscosity,  the  coefficient  of  heat  conduction,  the  specific 
heat  at  constant  pressure  and  the  enthalpy  per  unit  mass.  The  constants  7 ,  M  and  a  are  in 
turn  the  ratio  of  specific  heats,  the  Mach  number  and  the  Prandtl  number  defined  by 


_  f^oo  cpoo 

Jc 

•''OO 


where  K  is  a  gas  constant  and  =  ^7  is  the  sound  speed  in  the  free  stream.  Finally,  the 
function  a  in  the  equation  of  state  (2.9)  denotes  the  percentage  by  mass  of  the  mixture  which 
has  been  dissociated.  Later  in  §7  we  shall  give  the  expression  for  a  for  a  specific  dissociation 
model  used  in  our  discussion.  In  equations  (2.5)- (2.8),  the  operator  D/Dt  is  the  material 
derivative  and  it  has  the  usual  expression  appropriate  to  a  rectangular  coordinate  system. 

The  basic  state  is  given  by 

(«>  v>  w)  =  («(*. y),  Kx>y),  o),  t  =  f(x,y), 


p  =  p(x,y)>  p  =  p(x>y)-  (2-io) 

By  substituting  (2.10)  into  the  governing  equations  (2.4)-(2.9)  it  is  straightforward  to  obtain 
the  reduced  equations  satisfied  by  the  basic  state.  The  reader  is  referred  to  the  book  by 
Stewartson  (1964)  for  a  detailed  discussion  of  these  basic  state  equations.  If  we  define  the 
Howarth-Dorodnitsyn  variable  y  and  a  similarity  variable  77  by 

rv  y 

y=  pdy  and  77  =  -==,  (2.11) 

Jo  v2x 

then  the  continuity  equation  is  satisfied  if  u  and  v  are  written  as 

u  =  f\rj),  v  =  -L[-i/(77)  +  f'(v)  r  \drj\.  (2.12) 

V2x  p  Jo  P 

Here  the  functions  f(rj)  and  T(rj)  must  satisfy 


ff"  +  (m")'  =  0,  (2.13) 

\{pkf')'  +  bff  +  M(7  -  1  )M2p(f")2  =  0,  (2.14) 

if  the  x-momentum  and  energy  equations  are  to  be  satisfied.  These  equations  must  then  be 
solved  such  that  /,/'  vanish  at  the  wall,  f',T  =  1  at  infinity  and  either  T'  =  0  or  T  specified 
at  the  wall.  The  y-momentum  equation  gives 
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to  leading  order  so  that  p  =  p(x).  In  our  following  analysis,  we  assume  that  there  is  no  pressure 
gradient  along  the  stream  wise  direction  and  therefore  we  can  take  p  =  constant.  Equation  (2.9) 
then  reduces  to 

[1  +  a{T)}pT  =  1.  (2.15) 

Note  that  in  obtaining  (2.13)  and  (2.14)  we  have  not  made  any  constitutive  assumptions  so 
that  they  are  valid  for  dissociated  gases  (to  be  discussed  in  §7)  as  well  as  for  undissociated 
ideal  gases  (to  be  discussed  in  §3- §6). 


3  The  perturbation  equations 


We  first  assume  that  the  fluid  is  an  ideal  (one  component)  gas  undergoing  no  dissociation  so 
that  a  =  0.  Then  we  can  assume  that  (i).  the  specific  heats  are  constants;  (ii).  the  coefficient 
of  heat  conduction  is  linearly  related  to  the  shear  viscosity  and  (iii).  the  enthalpy  h  is  given 
by  h  —  cpT.  These  assumptions  lead  to  the  results 


k  =  p,  Cp  =  1,  p 


1 

f ' 


(3.1) 


(Note  that  all  of  these  quantities  have  been  non-dimensionalized).  Then  the  basic  equations 
(2.13)  and  (2.14)  simplify  to 

//"  +  (f/")'  =  0,  (3-2) 

-(£?')'  +  ff  +  (7  -  1  )M2hf"?  =  0.  (3.3) 

o  1  1 

These  two  equations  can  then  be  solved  if  we  make  an  constitutive  assumption  about  the 
viscosity  p.  In  the  previous  paper,  Hall  and  Fu  (1989),  we  used  Chapman’s  viscosity  law.  Here 
we  use  Sutherland’s  viscosity  law,  the  dimensionless  form  of  which  is  given  by 

f3/2 

p  =  (1  +  — r>  (3-4) 

where  m  is  a  constant.  Equation  (3.4)  is  exact  in  the  sense  that  it  is  derivable  from  the  kinetic 
theories  of  gases  (see  Chapman  and  Cowling  (1970)  for  a  discussion  of  its  validity,  also  compare 
(3.4)  with  (7.9)).  At  high  Mach  numbers  we  know  from  the  work  of  for  example  Freeman  and 
Lam  (1959)  that  the  basic  state  splits  up  into  two  distinct  regions.  Near  the  wall  a  boundary 
layer  forms  in  which  the  downstream  velocity  approaches  it’s  free-stream  value  of  unity  whilst 
the  temperature  decreases  from  it’s  value  at  the  wall  and  ultimately  decays  algebraically  at 
the  edge  of  the  layer.  In  the  next  region  this  algebraic  decay  is  taken  up  and  the  temperature 
then  approaches  exponentially  the  free-stream  value  of  unity. 

As  mentioned  above  an  explicit  analytical  solution  for  the  equations  (3.2)  and  (3.3)  is  not 
possible.  However,  an  asymptotic  analysis  in  the  large  Mach  number  limit  shows  that  the 
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boundary  layer  can  be  divided  into  two  regions:  an  inner  region  77  —  0(M  1//2)  and  an  outer 
region  77  =  0(1).  In  the  inner  region,  we  define  Y,T  and  F(Y)  by 

Y  =  M1/2  tj,  T  =  M~2T,  F{Y)  =  M1^2/.  (3.5) 

Then  equations  (3.2)  and  (3.3)  give 

(l  +  m)(-~  )'  +  FFtt  =  0,  (3.6) 

1  4-  fh  T'  -  ( F")2 

i-^(4 -)'  +  +  (7  -  1)(1  +  =  0  (3.7) 

O  s/t  vt 

to  leading  order.  These  equations  are  to  be  solved  numerically  subject  to  the  conditions 


F(  0)  =  F'(0)  =  0,  f  (00)  =  0,  F'(oc)  =  1, 


T'(0)  =  0  if  the  wall  is  thermally  insulated,  (3-8) 

T(0)  =  nTw  if  the  wall  is  under  cooling, 


where  Tw  is  the  wall  temperature  scaled  on  M2T(X  when  the  wall  is  thermally  insulated  and  n 
is  the  wall  cooling  coefficient. 

In  Fig.l  we  have  shown  the  results  of  our  numerical  integration  of  the  wall  layer  equations 
(3.6)  and  (3.7).  The  temperature  profiles  are  plotted  for  three  values  of  the  wall  cooling 
coefficients:  n=0.2,  0.6  and  1.0  and  were  calculated  with  7  =  1.4,  a  -  .72,  fh  =  .509.  The 
asymptotic  profile  for  large  Y  given  by  equation  (3.10)  is  also  plotted  there  for  comparison. 

For  large  Y,  equations  (3.6)  and  (3.7)  have  the  asymptotic  solutions 


F  =  Y  ~/3  + 


D 

(y  -  py/*  +  "  ‘  ’ 


(3.9) 


T  = 


9(1  +  fh)2 


{Y-py 


+ 


(3.10) 


where  both  (3  and  D  are  to  be  determined  by  a  numerical  calculation.  The  numerical  values 
of  (3  corresponding  to  four  values  of  the  wall  cooling  coefficients  are  listed  in  Table  1  together 
with  the  values  of  F"(0),  T'(0)  and  T( 0). 


Table  1 


n  =  0.2 

n  =  0.4 

n  =  0.6 

n  =  0.8 

n  =  1 

F"(  0) 

0.1517 

0.1997 

0.2317 

0.2560 

0.2758 

f'(  0) 

1.8192  x  10"2 

1.7976  x  10~2 

1.3909  x  10"2 

7.6872  x  lO"2 

0 

f(0) 

3.286  x  10-2 

6.572  x  10"2 

9.858  x  10“2 

1.3144  x  lO"1 

1.643  x  10 

0 

3.1808 

2.8366 

2.6301 

2.4840 

2.3721 
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The  asymptotic  expressions  (3.9)  and  (3.10)  imply  that  in  the  region  77  =  0(1), 

,  „  0  ,  hn) 

T  =  f(r,)+.... 

On  substituting  (3.11)  and  (3.12)  into  (3.2)  and  (3.3),  we  obtain  to  leading  order 


(1  +  rh) 


tTZ'"  +  ’""  =  0' 


1  +  m  /  Vf  „ 

- - T'  +  r?T'  =  0. 

o  \T  +  m  J 

These  two  equations  are  to  be  solved  numerically  subject  to  the  matching  conditions 


V~*0,  T 


9(1  +  rh)2  1 


and  the  conditions  at  infinity 


/'( 00)  =  0,  f(oo)  =  1. 


(3.11) 


(3.12) 


(3.13) 


(3.14) 


(3.15) 


(3.16) 


The  result  from  such  a  numerical  calculation  is  shown  in  Fig. 2.  In  this  figure  the  asymptotic 
result  is  the  one  given  by  (3.15). 

It  can  be  shown  from  (3.13)  and  (3.14)  that 


f  +  rh 

Vf 


(-T’,)1/<T, 


(3.17) 


so  that  after  (3.14)  has  been  solved  numerically,  the  function  f"{rj)  can  be  computed  easily 
from  this  equation.  Also,  we  note  that  whilst  the  solution  of  (3.14)  is  independent  of  the  inner 
region  solution  and  thus  of  the  conditions  at  the  wall,  the  function  /  is  dependent  on  the  inner 
region  solutions  through  the  matching  constant  D. 

We  now  assume,  as  in  Hall  and  Fu  (1989),  that  the  Sow  is  perturbed  to  spanwise  periodic 
stationary  vortex  structure  with  constant  wavenumber  a.  The  linearized  stability  equations 
for  these  Gortler  vortices  are  then  found  by  linearizing  the  Navier-Stokes  equations  (2.4)— (2.9) 
about  the  basic  state  and  retaining  the  leading  order  terms  in  the  high  Reynolds  number  limit. 
We  obtain 


y(«ff.  +  +  (Aa2  +  -  ( fiUv)v  +  luyV 

~  +  Vfiy)  +  (AUy)v}  T  ~  AW  =  °» 


(3.18) 


8 


i(vx  +  kuG)U  +  2-hUx  -  ~ P'Uxy  -  fixUy  +  i( uVx  +  vVy)  +  ( p.a 2  +  ^)V  -  £Fy)v 

+  Py  -  [j^{uVx  +  VVy  +  i/C Gu2)  +  ^ jluXy  -  ^flyVv  +  ^(A®v)v  +  A xUy]T 

2  4  2  1 

-  JlUyTx  -  [--fiux  +  -fLVy]Ty  +  -ftyiaW  -  -iajXWy  =  0,  (3.1?) 

1  _  _  1  2__ 

[LxiaU  +  -fiiaUx  +  [lyiaV  +  -fiiaVy  -  iaP  -  -p.(ux  +  vy)iaT 

<5  O  «) 

-  J,(uWx  +  vWy )  -  ^na2W  +  {pWy)y  =  0,  (3.20) 

jij(uTX  +  VTy)T  —  ^2  (ux  +  ®y)r  -  -^(uTx  +  VTy) 

+  f(Ux  +  Vy)  -  ±(TXU  +  fyV)  +  ia(Z)  =  o,  (3.21) 

±TXU  -  2(7  -  l)M2p.UyUy  +  ItyV  +  l{uTx  +  vTy)  +  i a2T 

-  {^(*fx  +  t;Tv)  +  (7  -  1)M2£^  +  i(/ifv)tf  J  T  -  \iityTy  -  i(/*Tv)y  =  0.  (3.22) 

Here  fi  =  dfi/dT,  whilst  (U,V,W),  P  and  T  denote  the  vortex  velocity  field,  pressure  and 
temperature,  respectively.  Equations  (3.18)-(3.22)  differ  from  equations  (2.11a-e)  given  in 
Hall  and  Fu  (1039)  only  in  that  the  bulk  viscosity  is  taken  to  be  zero  here;  that  assumption  is 
actually  implied  in  that  paper. 

It  was  shown  by  Hall  (1982)  that  in  the  incompressible  case  the  neutral  curve  for  small 
wavelength  vortices  has  G  ~  a4  and  that  the  vortices  are  confined  to  a  layer  of  depth  a-1/2 
where  the  flow  is  locally  most  unstable.  Hall  and  Malik  (1989)  extended  this  approach  to  the 
above  system  for  M  =  0(1)  and  wrote 

G  =  g0a 4  +  a3  +  • 


They  found  that  the  leading  order  growth  rate  a1 6*  has  6*  given  by 

,U2Ty  UUy 


6'  =  £  +  (-  - 
<7^2  T3 


oT2 


)9o- 


(3.23) 


In  the  neutral  case,  8*  =  0  and  (3.23)  then  determines  the  neutral  Gortler  number  g0  as  a 
function  of  17.  The  most  unstable  location  77*  is  where  go  has  its  minimum.  In  Hall  and  Fu 
(1989),  it  is  found  that  whei  1  Chapman’s  law  is  used,  77*  moves  away  from  the  wall  as  the  Mach 
number  increases.  It  is  also  found  that  the  basic  state  temperature  is  0(Af2)  over  most  of 
the  boundary  layer  and  decreases  rapidly  to  its  free- stream  value  over  a  logarithmically  thin 
adjustment  layer  sitting  at  the  edge  of  the  boundary  layer.  It  is  in  this  thin  layer  that  77* 
lies  and  hence  where  go  has  it’s  smallest  order  of  magnitude  in  the  large  Mach  number  limit. 
Thus  it  is  concluded  that  the  thin  temperature  adjustment  layer  is  most  susceptible  to  Gortler 
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vortices.  From  the  preceding  discussion  in  this  section  we  see  that  when  Sutherland’s  law  is 
used,  the  logarithmically  thin  temperature  adjustment  layer  corresponding  to  Chapman’s  law 
is  now  replaced  by  a  more  complex  temperature  adjustment  layer  of  0(1)  thickness.  If  we 
still  use  (3.23)  with  6*  =  0  to  calculate  the  orders  of  go  in  the  inner  layer  rj  =  0(M -1/2)  and 
the  outer  temperature  adjustment  layer  rj  =  0(1),  we  find  that  go  =  0(M 15/2)  in  the  former 
and  go  =  0(1)  in  the  latter.  Hence  once  again  the  temperature  adjustment  layer  is  most 
susceptible  to  Gortler  vortices  with  wavenumber  of  order  one  or  larger.  It  should  be  noted, 
however,  the  above  conclusion  is  based  upon  a  large  wavenumber  argument.  In  section  6  we 
shall  show  that  the  wall  layer  rj  =  is  actually  of  order  M3/2  thickness  in  terms  of 

the  physical  variable  y.  Thus  Gortler  vortices  with  wavelength  comparable  with  the  boundary 
layer  thickness  must  be  trapped  in  the  wall  layer  and  have  a  =  0(M-3/2).  It  will  be  shown 
in  section  6  that  this  wall  mode  has  neutral  Gortler  number  decreasing  monotonically  and 
has  the  centre  of  Gortler  vortex  activity  moving  towards  the  temperature  adjustment  layer  as 
the  wavenumber  increases.  Therefore,  the  minimum  Gortler  number  corresponds  to  the  mode 
trapped  in  the  temperature  adjustment  layer  and  the  latter  is  indeed  the  most  dangerous  mode 
when  the  whole  range  of  wavenumbers  are  considered.  It  should  also  be  noted  that  the  result 
go  =  0(1)  for  the  temperature  adjustment  layer  is  obtained  by  taking  the  large  Mach  number 
limit  of  the  0(1)  Mach  number  results.  By  doing  so  we  have  actually  missed  a  term  related 
to  the  curvature  of  the  basic  state  which  is  not  important  for  the  case  M  =  0(1)  and  a  >  1, 
but  is  important  in  the  large  Mach  number  limit.  As  we  shall  show  later  on,  the  curvature 
of  the  basic  state  produces  an  effective  Gortler  number  of  order  M 3/<2  in  the  absence  of  wall 
curvature  so  that  instability  can  not  occur  for  G  =  0(1). 

3.1  The  strongly  unstable  inviscid  mode 

Let  us  first  confine  our  attention  to  the  mode  trapped  in  the  temperature  adjustment  layer.  It  is 
easy  to  show  with  the  aid  of  expressions  (3.11)  and  (3.12)  that  in  this  temperature  adjustment 
layer, 

uv,  +  vv,  =  -p2_(BM3«  +  ^  -  ,2f']  +  o(l),  (3.24) 

where 

B  d=  lim  M-3/2  (°°  T{rj)dr}.  (3.25) 

M— too  Jo 

An  investigation  of  the  y-momentum  equation  (3.19)  shows  that  the  Gortler  number  must  be 
of  order  M3/2  in  order  to  enter  the  leading  order  analysis.  Thus  we  write 

i«(x)G  =  G*{x)M3/2,  (3.26) 

so  that  for  a  given  constant  Gortler  number  G  we  compute  G*(x)  using 

G'(x)  =  ^k(x)GM'3/2.  (3.27) 
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(3.28) 


For  convenience,  we  also  define  another  function  Q(x )  by 

(?(I)  =  (2^37?  ’ 

so  that 

uvx  4-  vvv  +  \k{x)Gu 2  =  (G*  -  Q)M3/2  +  o(M3/2).  (3.29) 

With  the  use  of  this  relation,  we  can  deduce  from  the  perturbation  equations  (3.18)-(3.22) 
that 

V  =  0{M3'4T),  W  =  0(M3/4T),  P  =  0(M3/2T),  U  =  OiM^T),  (3.30) 
and  that  for  fixed  77, 

where 

Mi  d=  l/M^  +  5. 

We  therefore  look  for  asymptotic  solutions  of  the  form 

T  =  exp  (m3'4  J* (3(x)dx\  {T0(x,v)  +  77)  +  •••}, 


U  =  Mf1  exp  [m3'4  JX  P{x)dx^  {UQ(x,V)  +  M~3^4Ui(x,t])  +  •••}, 

V  =  M3/4  exp  (m3/4  JX 0(x)dxj  {Vo(x,v)  +  M~3/4Vi( x,V)  +  •••},  (3.31) 

W  =  M3/4  exp  [m3>4  fp(x)dx)  {W0(x,v)  +  M"3/4Wi(i,t/)  +  •••}, 


P=  M3/2exp  (m3/4  jX  P{x)dxSj  {Po(x,v)  +  M~3/4Pi(x,  77)  +  ••■}, 

where  the  spatial  amplification  has  been  taken  care  of  using  the  WKB  method  and  /3(i)  is 
the  local  growth  rate  to  be  determined.  On  substituting  (3.31)  into  (3.18)— (3.22)  and  then 
equating  the  coefficients  of  like  powers  of  M,  we  obtain  a  hierarchy  of  equations.  To  leading 
order,  we  find  that  Vo  satisfies  the  differential  equation 


d*V0  iT'dVo  **  . 

WTW kTVo~7W’{  ~Q)TVo' 


whilst  To,  Wo  and  Pq  are  related  to  Vo  by 


„  _  .  ___  _  1  dV0  p  _  /V tedVp 

T°~  V2iT(3V°'  taW°~  \/2zT  dr)  ’  P°  T2k2  dr]  ’ 


(3.32) 


(3.33) 


and  Uo  does  not  appear  in  our  leading  order  analysis.  Here  k  d=‘  \/2®a  is  the  local  wavenumber. 
Equation  (3.32)  subject  to  Vq  vanishing  at  77  =  0, 00  is  a  Sturm- Liouville  problem  which  has 
solutions  if 

(G«  -  Q)f 

P  ~  ' 
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This  means  that 


/32  >0  if  G*  >Q,  /32  <0  if  G*  <  Q, 


since  T‘  <  0.  It  then  follows  that  neutral  stability  (/?  =  0)  occurs  at  the  position  x  —  xn  where 


G*(xn )  =  Q(xn) 


(3.34) 


at  zeroth  order.  Therefore,  in  view  of  the  definitions  (3.27)  and  (3.28),  the  neutral  Gortler 
number  has  the  expansion 


G  =  — — — - -tjtM3/2  +  higher  order  correction  terms. 

K(x„)(2xn)3/2 


(3.35) 


For  the  rest  of  this  paper  we  shall  take  Gn  to  be  the  zeroth  order  approximation  to  the  “critical 
inviscid  Gortler  number”,  thus  G n  is  obtained  by  retaining  only  the  first  term  on  the  right 
hand  side  of  (3.35).  An  important  point  concerning  (3.35)  is  that  the  first  term  on  the  right 

—  3 

hand  side  is  independent  of  xn  if  the  wall  curvature  varies  like  x  »  ;  in  the  latter  situation 
nonparallel  effects  dominate  and  the  vortex  growth  rate  is  smaller.  Thus  to  determine  the 
higher  order  correction  terms  to  the  neutral  Gortler  number,  we  have  to  distinguish  two  cases, 
namely,  (i).  k(x)  =  (2x)~3/2  and  (ii).  k(x)  ^  (2x)-3/2.  They  will  be  treated  separately  in  the 
next  two  sections. 

On  the  other  hand,  equation  (3.32)  can  also  be  interpreted  as  an  eigenvalue  problem  which 
determines  the  growth  rate  (3(x)  at  a  given  value  of  x  corresponding  to  any  wavenumber  k. 
The  appropriate  boundary  conditions  are  deduced  as  follows.  As  77  — ►  00,  T  — ♦  1  and  equation 
(3.32)  reduces  to 

d2V0  12  t/  n 
— —  -  k'VQ  =  0, 


so  that  Vo  ~  exp (-kr))  and  the  asymptotic  condition 

?p-  +  'kVo  =  0 

OTJ 


(3.36) 


should  be  imposed  at  “infinity”.  As  77  — ►  0,  T  — >  A2 jrf  where  from  (3.15b)  A  =  3(1  +  in)/ a. 
Equation  (3.32)  reduces  to 


d2V0  ,  80VO  k2A*  ( G'-Q)k 3  4 A 

- L - Vo  = - - 


dr/2  77  dr} 


\f2xfl2 


T 


which  has  the  solution 


Hence  the  asymptotic  condition 


T/  / 

V0~exp(--^-). 
dr]  r f 


(3.37) 


(3.38) 


(3.39) 


should  be  imposed  at  “zero”  (which  is  taken  to  be  some  small  value  in  numerical  calculation). 
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The  eigenvalue  problem  (3.32),  (3.36)  and  (3.39)  was  solved  numerically  by  employing  a 
fourth  order  Runge-Kutta  method.  In  Fig. 3,  we  have  shown  the  dependence  of  the  growth  rate 
on  the  local  wavenumber.  The  plot  clearly  shows  that  as  k  — »  0,  P2  — *  0  whilst  as  k  — *  oo, 
P2  — >  constant.  These  features  are  also  borne  out  by  the  asymptotic  analysis  given  in  the 
following  subsection.  The  inviscid  mode  we  have  described  above  therefore  has  growth  rate 
proportional  to  Mi  and  we  refer  to  it  as  the  strongly  unstable  inviscid  mode.  We  note  that 
when  G*  =  Q  the  growth  rate  vanishes.  In  this  case  it  is  necessary  to  look  for  evolution  of  the 
vortices  on  a  shorter  lengthscale  in  the  streamwise  direction;  that  problem  will  be  addressed 
later  in  this  paper  and  we  shall  refer  to  the  inviscid  mode  in  that  regime  as  the  near  neutral 
inviscid  mode. 


3.2  The  small  and  large  wavenumber  limits  of  the  strongly  unstable  inviscid 
mode 


First,  we  note  that  as  k  — ►  0,  we  are  approaching  the  scalings  for  the  wall  mode  (see  §6)  which 
has  wavenumber  k  ~  0(M~3/2)  and  which  is  trapped  in  the  wall  layer.  Thus  Gortler  vortices 
are  appropriately  governed  by  (3.37).  The  solution  (3.38)  shows  that  vortices  decay  to  zero  in 
the  thin  layer  r]  =  0(k1!3)  near  the  wall.  This  is  also  verified  by  the  obvious  shift  to  the  left  of 
the  first  mode  eigenfunctions  in  Fig.4  with  decreasing  k.  It  can  then  be  deduced  from  (3.37) 
that  P2  has  to  be  0{k )  in  order  to  enter  the  leading  order  analysis. 

Next,  in  the  large  k  limit,  a  WKB  analysis  of  (3.32)  shows  that  Gortler  vortices  will  be 
trapped  in  an  0(k~ 1/?2)  thin  layer  centred  at  Tf  =  rj*  where  P2  has  a  maximum.  Thus  we 
introduce  a  new  variable  £  by 

C^/^-r/*),  (3.40) 

expand  ( G *  -  Q)l(y/2xP2)  as 

=  X0  +  k~l'2X1  +  k-1  X2  +  •  •  • ,  (3.41) 

and  look  for  solutions  of  the  form 


V0(x,V)  =  V0°(x,0  +  k-l'2Vj(x,Q  +  k-'V2(x,0  +  •  •  • . 


(3.42) 


On  substituting  (3.40)-(3.42)  into  (3.32),  equating  the  coefficients  of  like  powers  of  k,  and 
then  solving  the  resulting  set  of  equations,  we  find  that  to  leading  order,  Aq  is  determined  as 


Aq 


fV) 

T'{v*Y 


(3.43) 


At  order  k  1/2,  Aj  is  determined  as  Ai  =  0  if  we  insist  that  P2  attains  its  maximum  at  tj  =  77*. 
At  order  A:-1,  Vj  is  found  to  satisfy  the  parabolic-cylinder  equation 


d2V$ 

d? 


-  4  W  -  o.Vg  =  0, 


(3.44) 
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where 


S=W){2?^)^(^}",/2-  P'45) 

If  we  impose  the  condition  that  the  disturbance  is  confined  to  the  k~ 1^2  layer  we  must  choose 
a  =  -(1/2)  —  3,  where  s  is  a  non- negative  integer.  This  condition  determines  the  infinite 
sequence  of  eigenvalues 


A2  =  A2 ,  =  Aq 


l  +  2s  /  d2  1 

-2T'(v*)  V  drj*3^X0  ' 


The  eigenfunctions  corresponding  to  the  eigenvalue  A2*  are  given  by 

^o°  =  V0°  =  e-tfHe.tf), 


(3.46) 


(3.47) 


where  He,(£)  is  a  Hermite  polynomial. 

With  the  aid  of  the  numerical  values  for  T  given  earlier  in  this  section,  we  find  that 

=  0.5786(1  -  0.529?!^  +  •••},  (3.48) 

and  that  77*  =  2.3228.  Both  these  numerical  values  and  the  concentration  of  vortices  in  an 
0(fc-1/2)  region  are  confirmed  by  the  numerical  results  shown  in  Fig.3  and  Fig. 5. 

In  closing  this  section,  we  note  that  results  given  here  for  hypersonic  flows  are  in  sharp 
contrast  with  related  results  for  incompressible  flows.  In  a  recent  paper,  Denier,  Hall  and 
Seddougui  (1990)  have  discussed  the  spectrum  of  the  large  Gortler  number  eigenvalue  prob¬ 
lem.  It  was  found  that  the  inviscid  Gortler  vortex  eigenvalue  problem  has  an  exact  solution 
with  the  spatial  growth  rate  increasing  monotonically  from  zero  and  tending  to  infinity  at 
large  wavenumbers.  In  the  high  wavenumber  limit  viscous  effects  become  important  when  the 
wavenumber  is  0(G® )  and  a  maximum  growth  rate  is  achieved  in  that  regime  with  the  growth 
rate  tending  to  zero  in  an  0(G«)  regime  as  discussed  by  Hall  (1982).  At  small  wavenum¬ 
bers  the  growth  rate  tends  to  zero  and  the  vortices  spread  out  above  the  boundary  layer.  In 
fact  when  the  wavenumber  is  0[G~T)  an  eigenvalue  problem  related  to  that  for  Tollmien- 
Schlichting  waves  is  recovered  but  at  even  smaller  wavenumbers  nonparallel  effects  dominate 
and  the  problem  must  be  solved  numerically  as  in  Hall  (1983).  The  major  difference  we  have 
found  above  for  hypersonic  Gortler  vortices  is  that  the  growth  rate  tends  to  a  finite  value  as 
the  local  wavenumber  tends  to  infinity;  we  shall  see  later  that  this  has  a  significant  effect  on 
the  way  in  which  viscous  effects  come  into  play  at  high  wavenumbers. 


4  Neutral  instability  with  «(x)  ~  (2x)-3/2 

We  now  proceed  with  the  determination  of  the  higher  order  correction  terms  in  the  neutral 
Gortler  number  expansion  (3.35).  In  the  case  when  the  curvature  k(x)  =  (2®)-3/3,  G*(x )  = 


<J(i)  and  the  0(M3/2)  term  on  the  right  hand  side  of  (3.29)  vanishes  for  all  x.  The  implication 
of  this  is  that  for  this  special  distribution  the  curvature  of  the  basic  state  is  exactly  counteracted 
by  wall  curvature  over  an  0(1)  interval  in  x,  in  the  more  general  case  that  is  only  the  case  over 
an  asymptotically  small  interval.  An  investigation  of  the  perturbation  equations  (3.18)-(3.22) 
with  the  aid  of  the  equations  (A1)-(A5)  given  in  the  Appendix  then  reveals  that  the  neutral 
Gortler  number  must  expand  as 

G  =  2BM3/2  +  G  +  o(l),  (4.1) 


and  that  the  perturbation  quantities  have  relative  orders 

u  =  °(j^H  T=0(V),  W  =  0(V),  P  =  0(V), 

where  G  =  0(M°)  is  to  be  determined.  We  therefore  look  for  the  following  form  for  the 
solutions  for  (3.18)-(3.22): 

U  =  W1&{x'v)  +  '"i  v  =  V(x,v)  +  --->  W  =  W(x,r,)  +  ...,  (4.2) 

p=  ~ZrHx>  7/)  +  ”->  T  =  V2xO(x,r])  +  -.-, 

V2x 

where  the  insertion  of  the  factor  \f2x  is  purely  for  convenience. 

On  substituting  (4.2)  into  (3.18)-(3.22)  and  then  equating  the  coefficients  of  like  powers 
of  M,  we  obtain,  to  leading  order, 


£0#  dU_ 
dryT  07/ j  + 77  07? 


=  V2x 


\TV  +  [ff  077 (  f  "  f  077 


d±  =  OLv  +  77^  -  aVf V  + 

077  T  +VdV  M  +  3  077lr  077  ; 

+irr,2T'  _  EiL  4.  ir  _  4t  9 

+f[v  +  2G  3Tdrj(  f  )J 

2.  9W  4  „  SW 

-  -infr  ike  -v-g^  +  jM  TW  +  2*-^-, 

0  ,V,  r)d§  ,T-  „  r,T'J  ,  2 xd§ 

077(7)+  y 077  -  1  W  +  (  +  T  )f  +  f  dx' 


(4.3) 

(4.4) 

(4.5) 

(4.6) 
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(4-7) 


Here  k  d=  V2xa.  It  can  be  seen  that  (4.4)-(4.7)  are  independent  of  U  and  the  latter  is 
determined  from  (4.3)  after  (V,  W,  f ,  P)  have  been  determined.  Since  these  leading  order 
perturbation  equations  are  parabolic  with  respect  to  the  variable  x,  they  have  to  be  solved 
by  specifying  the  perturbation  quantities  at  a  given  upstream  position  and  then  marching 
downstream.  We  therefore  expect  that  neutral  stability  would  depend  crucially  on  what  initial 
conditions  we  impose  and  where  we  impose  them.  However,  before  we  present  our  numerical 
solutions  of  these  equations,  we  first  consider  a  special  case:  the  large  wavenumber  limit  for 
which  a  simple  asymptotic  solution  is  possible. 


4.1  Large  wavenumber  limit 

In  the  large  wavenumber  limit,  the  length  scale  over  which  vortices  vary  is  small  compared 
with  the  lengthscale  over  which  the  boundary  layer  grows.  Then  we  expect  that  nonparallel 
effects  do  not  come  into  our  leading  order  analysis.  This  is  indeed  the  case,  as  we  show  below. 

For  large  jfc,  vortices  are  confined  to  a  thin  layer  of  0(e1//2)  thickness  centred  on  tj  =  rj 
where, 

e  d=’  l/k,  (4-8) 

and  where  rf  is  the  most  unstable  position  to  be  determined  in  the  course  of  our  calculation. 
We  therefore  define  a  new  variable  (j>  by 

<f>  =  e~1/2(y 7-?7*). 

An  investigation  of  equations  (4.3)— (4.7)  shows  that  when  c  — *  0, 

0  =  0(e2V),  W  =  G(e1/2F),  P  =  0{e~l‘2V), 

6  =  0(  l/<4),  ^  =  0(l/e!).  (4.9) 

Hence  we  look  for  the  following  form  of  asymptotic  solutions  for  (4.3)-(4.7) 

G  =  4(Go  +  «1A'<  ■  i  +  +  «3/2^3  +  •••). 

V  =  (Vo  +  e^Vi  +  •  •  -)E,  W  =  (61/2Wo  +  eWx  +  •  •  )E,  (4.10) 

P  =  (€-1/2P0  +  Pi  +  •  •  •),  9  =  (e20o  +  eS/29i  +  •  •  -)E , 

where  . 

E  =  ex pj^  J  +  f1/3A(^)  +  1 '  •)<#]  i  (4-11) 
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and  where  Vo,  V\  etc.  are  functions  of  <j)  and  x.  Note  that  E  here  represents  the  fast  variation 
of  the  perturbation  quantities  along  the  streamwise  direction  whilst  the  dependence  of  Vo,  Wo 
etc.  on  x  represent  the  slow  variation  of  perturbation  quantities  due  to  the  nonparallel  effect 
of  the  boundary  layer  growth.  In  effect  we  have  described  the  fast  variation  of  the  disturbance 
by  a  WKB  type  of  expansion  in  the  streamwise  direction.  Here  we  are  only  concerned  with 

neutral  stability,  so  we  set  Qo  =  Pi  =  02  =  0.  On  substituting  (4.10)  into  (4.3)-(4.7)  and  then 

equating  the  coefficients  of  like  powers  of  e,  we  obtain  a  hierarchy  of  matrix  equations.  To 
leading  order,  we  have 

a  _  _  1  dV0 

e°--Kr!v°‘ 

_  -  fiodVo  _  2/lofo 

0  To  d<f>  ’  G°  oT\  '  ^4'12^ 

where  To  =  T(rj*),Ti  —  T'( rf)  and  /Zo  =  £(To).  To  next  order,  we  deduce  that 

G\  =  0.  (4.13) 

At  next  order  we  find  that  Vo  must  satisfy  the  parabolic-cylinder  equation 

m  =  o.  (4.i4) 

Here 

C  =  ^,  (4.15) 

If  we  impose  the  condition  that  the  disturbance  is  confined  to  the  e1/2  layer  we  must  choose 


d  =  -1/2  -  s,  s  =  0, 1,-  •  • , 

The  smallest  G2  corresponds  to  3  =  0  and  we  then  have  from  (4.15b)  that 

^  _  (3God2Go\1/2 

G’  =  UiFlv  J. ' 


(4.16) 


(4.17) 


The  centre  of  vortex  activity  r f  is  determined  by  the  condition  that  Go  attains  its  minimum 
there: 


<isf w  -  °- 


(4.18) 


After  solving  (3.14)  numerically  for  the  basic  state  temperature  T,  we  then  use  (4.12d)  and 
(4.18)  to  determine  rf ,  and  (4.12d)  and  (4.17)  to  determine  Gq  and  G 2.  We  find  that 


rf  =  3.001,  Go  =  15.4834,  G2  =  34.3175, 


so  that 


G  =  15.4834&4  +  34.3175P  +  •  •  • .  (4.19) 

Finally,  we  remark  that  the  above  analysis  is  valid  as  long  as  the  local  wavenumber  k  =  \/2xa 
is  large.  This  means  that  the  far  downstream  evolution  of  Gortler  vortices  can  always  be 
described  by  the  above  theory. 
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4.2  Numerical  results  for  k  0(1) 

When  the  wavenumber  is  0(1),  the  perturbation  equations  (4.3)-(4.7)  are  parabolic  with 
respect  to  x  and  they  have  to  be  solved  numerically  subject  to  some  initial  conditions  imposed 
at  some  upstream  position  x.  For  computational  purpose  it  is  convenient  to  eliminate  P  and 
W  among  (4.4)-(4.5).  After  some  manipulation,  we  obtain 


where 


ldAV  2x  cPV  4 xf'd2V  k2T2  n  dV  A 

— - 1 - - - 2x —  =  >  a,. 

T  drf  p  dxdr,2  pT  dxdr)  \i  dx  — 


cT  dr?  ~  1  oT*  cT  11  dr,  T  Sx 

a,  = 

a,  =  P^vfF  -  vf  + 

P  n  dr, 

a3  =  -k2  /(I  _  I)f '  +  rff"  -  ^  -  ^-(p' 7}f'  +  pf  +  prjf") 

[  CT  1  H 

+^  +  HL  +  ^HL  +  La\e, 

P  p  p  2y.  ) 

a*  =  -  £)  -  **<  J  +  ±)(T"T  -  1"*)] 

-Sa(|  -  ^)f'!  +  iF(|f »  +  T"-H2) 


+  ^-r,Tf 

P 


,-ifc4f3|(-F), 


3T'd3V  d  3?'  d2V  d2  T'  dV  d3 W3 

°5  “  "  T2  dr] 3  dry  T2  '  dr]2  dr,2^f2’  dr,  dr,3 


T2*2  f  2A^/  fr,  V  T'  1  ,  2xT 

\  p  P  T  k2pf 2  J  P 


(4.20) 


(4.21) 


d2Wi  n  f'  dWi 

1. 1,.. .  _ 2,x — 

dxdr,  p  dx 


k2T2n  dd 

0s  =  ___2x,_, 


U7  =  -2x  Tk2 


drjdx  ’ 


f  3T'  _  tfT  _  2A^,\  PW 

IT"  p  p  )d  v2’ 
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ag 


■( 


2 T2P  +  T  .1“.  +  (£  _  -  T")  -  ^T'l 

yfj,  fJK  T  ji  \  drj 


t]T‘  +  2f 


In  these  equations,  W3  and  W  are  given  by 

AT 


W3  = 


A  d20  +  ^2pT' 


ctT2  a??2 


aT2  aT3 


P^_sdO_  f_l 

ctT3^  \a^ 


Z)Ji2 

af  T  j  ct 


0, 


T?r  1  «97  TXr 
W  =  -f~^  +  W3- 


(4.22) 


(4.23) 


The  finite  difference  scheme  used  here  is  similar  to  that  used  by  Hall  (1983),  Wadey  (1990)  for 
incompressible  flows.  The  reader  is  referred  to  the  latter  papers  for  a  more  detailed  discussion. 
Our  implementation  of  the  numerical  scheme  is  as  follows.  We  first  specify  V",  U  and  0  at  a 
given  upstream  position  x.  Then  v;e  use  the  finite  difference  scheme  to  march  downstream  and 
thus  calculate  the  evolution  of  the  initial  perturbation  with  respect  to  the  streamwise  variable 
x.  The  position  of  neutral  stability  is  defined  as  the  place  where  a  certain  energy  measure  has 
zero  growth  rate.  We  use  the  following  five  energy  measures  to  monitor  the  energy  variation: 


E 1 =  f 

Jo 


°°  U2dy  =  V2x  r  U2Tdr], 
Jo 


E2  =  [ °° {V 2  +  W2)dy  =  /°°(F2  +  W2)fdv, 

Jo  Jo 

E3  =  /  \PIx0dy  —  2x  f  OfdTj,  £4  =  Umax, 
Jo  Jo 

5~Jo  T  dr,  dri' 


(4.24) 

(4.25) 

(4.26) 

(4.27) 


We  give  results  for  several  measures  of  the  vortex  strength  in  order  to  show  the  relatively  large 
variation  of  growth  rates  associated  with  different  flow  quantities.  Depending  on  the  initial 
conditions  and  on  which  energy  measure  we  use,  the  disturbance  can  either  grow  or  decay 
initially.  In  the  former  case,  there  is  only  one  neutral  position  corresponding  to  each  pair 
(a,  (5):  the  energy  will  reach  a  maximum  at  some  downstream  location  and  then  decay  to  zero 
monotonically.  The  corresponding  neutral  curves  only  have  right  branches.  The  region  on  the 
left  of  a  neutral  curve  is  unstable  and  the  region  on  the  right  is  stable.  In  the  latter  case,  if 
the  Gortler  number  is  large  enough,  there  are  always  two  neutral  positions:  the  energy  will 
reach  a  minimum  at  some  downstream  location,  then  grow  to  reach  a  maximum  and  finally 
decay  to  zero  monotonically.  The  corresponding  neutral  curves  have  distinct  right  and  left 
branches.  The  region  above  a  neutral  curve  is  unstable  and  other  regions  are  stable.  For  the 
special  curvature  case  discussed  in  the  present  section,  the  flow  is  always  stable  for  i>1,  since 
according  to  the  asymptotic  result  obtained  in  the  previous  subsection,  the  flow  is  neutrally 
stable  where  G  ~  k*  for  k  »  1  so  that  k(x)  must  increase  at  least  as  quickly  as  x1^2  if  the 
vortex  is  to  be  unstable  for  x  »  1. 
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At  a  neutral  position,  we  calculate  the  local  wavenumber  ax  and  local  Gortler  number  Gx 
defined  by 

ax  =  k  =  y/2xa,  Gx  —  /c(x)(2x)3,/2<j  =  G.  (4.28) 

By  fixing  a  and  varying  (5,  we  obtain  a  series  of  such  points  (aX)  Gx)  and  thus  plot  a 
neutral  curve.  In  our  calculation,  a  is  fixed  at  0.1;  the  step  lengths  along  x  and  rj  directions 
are  taken  to  be  0.1  and  0.05  respectively;  and  the  lower  boundary  of  the  0(1)  temperature 
adjustment  layer  is  taken  to  be  0.8  whilst  infinity  is  approximated  by  20.8. 

In  Fig.6-Fig.12,  we  have  given  the  results  of  our  numerical  calculations.  Fig.6  shows  five 
different  neutral  curves  obtained  when  we  use  the  five  different  energy  measures  (4.24)— (4.27) 
to  monitor  the  energy  growth.  Fig. 7  shows  three  different  neutral  curves  which  are  obtained 
when  we  impose  three  different  initial  conditions  at  the  upstream  position  Xq  =  20,  whilst  Fig.8 
shows  three  different  neutral  curves  which  are  obtained  when  we  impose  an  initial  condition 
at  three  different  upstream  locations.  In  two  of  these  graphs,  we  have  also  plotted  the  two 
term  asymptotic  result  (4.19).  As  we  expect,  although  these  neutral  curves  have  distinct  left 
branches,  their  right  branches  all  converge  to  the  unique  large  local  wavenumber  limit.  We 
can  see  from  Fig.8  that  as  the  initial  location  of  the  disturbance  moves  towards  the  leading 
edge,  the  neutral  curves  move  progressively  up  and  across  to  the  right.  This  is  contrary  to 
the  corresponding  results  found  by  Hall  (1983)  for  incompressible  flows.  Finally,  Fig.9-Fig.12 
show  the  typical  profiles  of  the  four  perturbation  quantities  as  they  evolve  downstream,  with 
the  initial  conditions  given  by  (iii)  in  Fig.6  and  Xo  =  20,(5  =  1000.  In  order  to  see  how  the 
the  centre  of  vortex  activity  evolves  downstream,  we  have  normalized  each  of  the  perturbation 
quantities  by  its  maximum.  It  is  clear  that  as  the  vortices  evolve  downstream  (and  thus  as  the 
local  wavenumber  increases),  they  become  more  and  more  concentrated,  which  agrees  with  the 
asymptotic  result  found  in  the  previous  subsection  that  in  the  large  local  wavenumber  limit, 
Gortler  vortices  are  trapped  in  a  thin  layer  of  depth  0(k~ 1/2)  centred  at  rj  =  3.001. 

5  Neutral  instability  with  k(x)  ^  (2x)-3/2 

When  the  wall  curvature  is  not  proportional  to  (2i)-3/2,  the  0(M 3/2)  term  on  the  right  hand 
side  of  (3.29)  only  vanishes  at  the  leading  order  neutral  position  and  its  effect  will  persist  in 
the  downstream  development  of  Gortler  vortices.  An  important  consequence  of  such  an  effect 
is  that  non-parallel  effects  will  be  important  over  a  larger  range  of  wavenumbers  than  was 
the  case  for  the  special  curvature  case.  Suppose  we  measure  the  order  of  the  wavenumber  by 
writing  it  as  a  =  0(Ma).  Then  we  will  show  in  this  section  that  non-parallel  effects  continue  to 
be  dominant  for  a  up  to  and  including  1/4.  For  a  >  1/4,  non-parallel  effects  become  negligible 
compared  with  viscous  effects  and  an  analytical  expression  can  be  obtained  for  the  second 
order  correction  to  the  Gortler  number  expansion.  This  second  order  correction  becomes  of 
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the  same  order  as  the  leading  order  term  due  to  the  curvature  of  the  basic  state  when  a  =  3/8. 
In  this  case  a  higher  order  correction  term  can  also  be  obtained. 


5.1  0(1)  wavenumber  regime — the  near  neutral  inviscid  mode 


In  the  0(1)  wavenumber  regime,  it  is  convenient  to  determine  the  stability  properties  by 
considering  the  evolution  of  Gortler  vortices  in  the  neighbourhood  of  the  leading  order  neutral 
position  xn  given  by  (3.34).  Thus  we  shall  fix  the  Gortler  number  as  given  by 


G  = 


_ — _ M3'2 

(2in)3/2  K(xn) 


(5.1) 


and  determine  the  second  order  correction  say  xn  to  the  neutral  position  xn  so  that  Gortler 
vortices  with  G  given  by  (5.1)  are  neutrally  stable  at  location  z„  +  xn.  Replacing  in  by  xn  -  xn 
in  (5.1)  then  gives  the  appropriate  expansion  of  the  Gortler  number  for  vortices  neutrally  stable 
at  x  —  3/fi  • 

It  can  be  shown  that  in  the  neighbourhood  of  xn,  the  second  term  in  the  expansion  of 
k(x)G/2  will  force  a  growth  rate  of  order  M1^2.  Hence  we  shall  consider  the  evolution  of 
Gortler  vortices  in  an  0(Af-1/2)  neighbourhood  of  x„  by  defining  a  new  variable  X  by 


X  =  (x-  x^M1'2, 


(5.2) 


and  look  for  asymptotic  solutions  of  the  form 

T  =  T0(X,  V)  +  ■  •  • ,  V  =  M^2Vo(X,  7?)  +  • . 

W  =  M^2W0(X,t,)  +  P  =  MPq(X,  77)  + 
Equation  (3.29)  becomes 


where 


uvx  +  vvv  +  -k(x)Gu2  =  EXM  +  o(M), 


p  def .  d{G'-Q) 
dx 


(5.3) 


(5.4) 


(5.5) 


Note  that  it  is  this  term  that  gives  rise  to  a  local  growth  rate  of  order  M1^2.  On  substituting 
(5.3ar-d)  into  the  perturbation  equations  (3.18)-(3.22),  dropping  higher  order  terms  and  then 
eliminating  Vo,  Wo  and  Pq  from  the  resulting  equations  in  favour  of  To,  we  find  that  Vq,Wq 
and  Pq  are  related  to  To  by 


V0  =  - 


ikW0  =  - 


y/2 xnT  dTo 
f  dX’ 
1  dV0 


Po=- 


T  dv  ’ 
y/2zn  d2Vp 
k2T2  dXdrj 2  ’ 


(5.6) 

(5.7) 

(5.8) 
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and  that  To  satisfies  the  equation 


d 2  [a2T0 

dx 2  \  dv2 


2 T"  dTo 
T'  dr] 


-[k2T2-TT\^ 


rptf 

ififn )  l-^0 


-lb2- 


EXT' 

V2x^ 


To  =  0, 


(5.9) 


where  k  d='  y/2x na.  We  can  interpret  (5.9)  as  the  turning  point  equation  associated  with  the 
breakdown  of  the  WKB  structure  in  x  of  the  expansions  (3.31),  indeed  the  evolution  equation 
(3.32)  is  retrieved  from  (5.9)  by  taking  X  to  be  large.  Since  this  breakdown  is  associated 
with  a  simple  zero  in  i  of  the  right  hand  side  of  (3.32)  we  expect  that  the  local  behaviour  of 
the  disturbance  should  now  be  expressed  in  terms  of  Airy  functions.  Equation  (5.9)  admits 
separable  solutions  of  the  form 

T0{X,r,)  =  ^(X)rP(rj),  (5.10) 


with  0  and  Tp  satisfying 

j>"(X)-uX4>(X)  =  0,  (5.11) 

1  rpH  P'V1 

nv)  -  -jt-av)  - 1 k2f 2  -  Tf\±  -  fosmv)  -  *2  •  =  °»  (5-12) 

where  the  separation  constant  w  is  to  be  determined  by  solving  the  eigenvalue  problem  (5.12) 
subject  to  appropriate  boundary  conditions.  By  a  simple  substitution  z  =  Xu;1/3,  equation 
(5.11)  reduces  to  the  standard  form  of  Airy’s  equation  W"{z)  -  zW(z)  =  0  which  has  two 
independent  solutions  Ai(z)  and  Bi(z),  so  the  solution  of  (5.11)  is  given  by 


4>{X)  =  aAi{u>1/3X)  +  65i(u>1/3X), 


(5.13) 


where  a  and  b  are  two  constants  to  be  determined  by  initial  conditions. 

To  solve  (5.12),  we  first  note  that  in  the  large  wavenumber  limit,  equation  (5.12)  takes 
the  same  form  as  equation  (3.32).  Therefore,  the  solution  of  (5.12)  can  be  written  in  terms  of 
Hermite  polynomials  as 

Tp{v)  =  V’*  =  e~«€S  •  He.(Z)  (5.14) 


and  from  (3.41)  the  eigenvalue  u;  expands  as 


(5.15) 


where  \q,  and  £  are  defined  in  turn  by  (3.43),  (3.4b)  and  (3.45a). 

In  the  0(1)  wavenumber  regime,  equation  (5.12)  has  to  be  solved  by  a  numerical  integration, 
and  in  general  an  infinite  number  of  eigenvalues  uj,(s  =  0, 1,  •  •  •)  and  eigenfunctions  ip,  can  be 
obtained.  Then  the  general  solution  of  (5.9)  can  be  written  as 


OO 

To  =  {a.Ai(u.X)  +  btBi{ujtX)}y>.(ri), 


i=0 


(5.16) 
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where  a,  and  b 4  are  constants  to  be  fixed  by  initial  conditions  at  X  =  0.  From  (5.6),  Vo  is 

given  by 

Vo  =  {aaAi'(u>,X)  +  btBi'{u.X)}  (5-17) 

*  s= 0 

It  is  clear  that  once  To(X,ri)  and  Vo(X,tj)  are  specified  at  X  =  0,  the  coefficients  (a4,64)  and 
hence  the  evolutionary  behaviour  of  the  perturbation  field  (Vo,  Wo,  To,  Pq)  will  be  completely 
determined. 

The  correction  term  to  the  neutral  position  can  be  defined  as  the  position  where  a  certain 
energy  measure  has  zero  growth  rate.  It  is  obvious  that  such  a  position  would  depend  upon 
what  initial  conditions  we  impose  at  X  =  0  and  what  energy  measure  is  employed  to  monitor 
the  energy  growth.  In  principle  then  it  is  an  easy  matter  to  determine  the  local  neutral  position 
associated  with  any  initial  perturbation,  we  note  however  that  before  growth  of  the  vortices 
occurs  they  will  have  an  oscillatory  behaviour  in  X  since  both  Airy  functions  are  oscillatory 
on  the  negative  real  axis.  Clearly  this  occurs  because  the  boundary  between  instability  and 
stability  is  controlled  by  inviscid  effects  in  this  regime,  there  is  no  counterpart  to  this  result  in 
the  behaviour  of  Gortler  vortices  or  for  that  matter  Tollmien-Schlichting  waves  in  incompress¬ 
ible  flows.  We  further  note  that  appropriate  forms  for  the  initial  conditions  can  be  obtained 
from  the  receptivity  problems  associated  with  wall  roughness  or  free  stream  disturbances,  see 
Denier,  Hall,  and  Seddougui  (1990)  and  Hall(1990).  We  merely  note  in  passing  here  that  it  is 
reasonable  to  expect  that  the  type  of  mode  discussed  above  is  more  likely  to  be  stimulated  by 
free-stream  disturbances  since  the  effect  of  wall  roughness  is  diminished  by  the  wall  layer  over 
which  the  wall  roughness  must  diffuse  before  reaching  the  unstable  adjustment  layer. 

This  is  certainly  typical  of  the  evolution  of  Gortler  vortices  in  growing  boundary  layers.  In 
the  present  problem,  non-parallel  effects  dominate  in  the  evolution  of  Gortler  vortices  mainly 
through  the  0(Af3/2)  curvature  of  the  basic  state.  As  we  increase  the  wavenumber,  viscous 
effects  will  gradually  come  into  play  in  the  evolution  of  Gortler  vortices  and  nonparallel  effects 
will  become  less  important.  In  the  following  subsection  we  consider  wavenumbers  of  order 
M1/4.  This  is  the  maximum  order  at  which  nonparaliel  effects  are  dominant.  We  shall  show 
that  when  the  wavenumber  is  increased  further  above  this  order,  nonparallel  effects  become 
negligible. 

5.2  The  0(M1/4)  wavenumber  regime — the  nonparallel  viscous  mode 

When  the  wavenumber  reaches  the  order  M1/4,  the  streamwise  lengthscale  is  still  0(M-1/2) 
(implying  that  the  local  growth  rate  is  0(M 1^2)),  but  the  vertical  lengthscale  becomes  of 
0(M~1/8)  (as  we  expect  that  vortices  would  be  trapped  in  an  0(M-1/8)  thin  layer).  We 
therefore  define  two  new  independent  variables  X  and  £  by 

X  —  (x  -  xn)M1,/2,  £  =  (,?-  (5.18) 
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where  77*  is  the  center  of  vortex  activity.  For  convenience,  we  define  a  small  parameter  e  by 

(5.19) 

An  analysis  of  the  perturbation  equations  (3.18)-(3.22)  shows  that 

V  =  0(e~AT),  W  =  0(e~3T),  P  =  0{e~hT).  (5.20) 

We  now  assume  that  T  =  0(1)  and  scale  (V,W,P)  by  (e-4,  e-3,  e-5).  Then  by  neglecting 
terms  of  relative  order  higher  than  c2  in  the  perturbation  equations  (3.18)-(3.22),  we  obtain 

dV  EX_  ,  4  d  ,udV. 


+  h2,TV-?!T-S  4  d(tdV) 

+  M  TT  3  •  (2xn)  f 

+  £2._L.9L_e2^_dW=  o 

v^2zn  dt,  3v/2xn  d£  ’ 


(5.21) 


1  dW  -  4  -  -  p,  dV 

¥w-TP+rTW+3jki(+‘ 


2  1  B  p  9* 

(2xn)a2dCf  d(  '  ’ 

(5.22) 

1  dV  1  dT  T'  Tr 

(5.23) 

w  = 

V2^Td(+e  TdX+e  V^f2V' 

V 

T  dT  a2 /if2  2  f  d  pdT 

V  dX  aV  2+6  (2  xn)oT'  d(  d( 

(5.24) 

where  W  d='  iaW,  a  d=‘  M_1/4a  and  where  we  have  used  the  same  notation  to  denote  the  scaled 
perturbation  quantities.  Elimination  of  V,  W,  and  P  in  favour  of  T  among  these  equations  then 
gives 


d2T  „  U.2  *  9T  1.4  EXTi  _ 

+  (*  +  ~)a  +  -«  (5)  To  ?+  1  ~ — -  2 ^ 

0A2  47  dX  a  V%x^T0 


(5.25) 


where  f0  =  T(rj*),fi  =  T'(r]*),fio  =  /i(fo)  and  where  77*  is  chosen  at  higher  order  such  that  the 
vertical  structure  of  T  can  be  expressed  in  terms  of  parabolic  cylinder  functions  which  vanish 
at  77  =  ±00.  The  solution  of  the  above  equation  is  easily  found  to  be  expressible  as  the  product 

-X(l+l)SaHbT0 

of  an  Airy  function  multiplied  bye  j  .  As  before  the  Airy  function  grows  or  decays 
exponentially  for  large  X  and  is  oscilatory  on  the  negative  real  axis.  However  the  presence  of 
the  exponential  factor  now  leads  to  a  crucial  change  in  the  nature  of  the  streamwise  evolution 
of  the  disturbance.  We  refer  to  the  fact  that  the  exponential  factor,  induced  by  viscous  effects, 
now  means  that  in  the  stable  regime  the  disturbance  decays  exponentially  rather  than  oscilating 
as  was  the  case  previously.  This  result  is  consistent  with  the  usual  result  of  stability  theory 
that  inviscid  disturbances  change  from  being  oscilatory  to  being  exponential  in  character  when 
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instability  occurs;  viscous  instabilities  on  the  other  hand  are  exponential  in  nature  either  side 
of  the  stability  boundary.  The  local  neutral  position  can  only  be  obtained  by  specifying  an 
initial  disturbance  and  finding  where  the  solution  of  the  evolution  equation  begins  to  grow. 

If  we  let  d<l,  then  to  leading  order  equation  (5.25)  reduces  to 


d2T 
dX 2  + 


EXTj 

v/S\2 


T  =  0 


(5.26) 


which  matches  with  the  asymptotic  form  of  (5.9)  when  k  >  1. 

On  the  other  hand,  if  we  let  a  »  1,  then  to  leading  order  (5.25)  reduces  to 


EXT j 


+  -a4p.2T i2 

<y  i 


T  =  0, 


(5.27) 


which  shows  that  the  second  order  correction  term  to  the  leading  order  neutral  position  becomes 
independent  of  nonparallel  effects  and  is  given  by 


X  =  Xn  =  -^tiT*a\  (5.28) 

On  replacing  xn  by  xn  —  M~1/2Xn  in  (5.1)  and  expanding  the  resulting  expression  in  the 
neighbourhood  of  xn,  we  obtain  the  following  expression  for  the  Gortler  number  for  Gortler 
vortices  neutrally  stable  at  location  x  —  xn: 


G  = 


2BM3/2  4_ 

«(i„)(2In)V2+“®  + 


(5.29) 


where  _ 

-  _  2\/2®n/io?o 

9o  - - if, - 

ck0Ti 

The  expression  (5.29)  is  valid  for  wavenumbers  of  order  Ma,  1/4  <  a  <  3/8.  When  the 
wavenumber  reaches  the  order  M3^6,  the  second  order  correction  term  in  (5.29)  is  as  large  as 
the  first  term  and  a  more  accurate  asymptotic  expression  can  be  obtained.  That  situation  will 
be  discussed  in  the  following  subsection.  We  have  seen  above  that  in  the  0(M ? )  wavenumber 
regime  viscous  effects  come  into  play  and  modify  the  evolution  of  the  near  neutral  inviscid 
mode,  for  that  reason  we  refer  to  the  mode  in  this  case  as  the  nonparallel  viscous  mode.  We 
further  note  that  the  strongly  unstable  inviscid  mode  connects  directly  with  the  parallel  viscous 
mode  discussed  below.  Thus  we  have  shown  above  that  the  nonparallel  viscous  mode  connects 
with  the  near  neutral  inviscid  mode  in  the  vicinity  of  the  right  hand  branch  of  the  neutral 
curve.  Hence  the  initial  stages  in  the  evolution  of  the  right  hand  branch  of  the  neutral  curve 
are  governed  by  an  interplay  between  viscous  and  inviscid  effects.  We  further  note  that,  in 
view  of  the  limiting  form  (5.29)  valid  for  large  a,  as  the  wavenumber  increases  the  neutral 
Gortler  number  will  also  increase. 
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5.3  The  O(MH)  wavenumber  regime — the  parallel  viscous  mode 

When  the  wavenumber  becomes  of  order  it/3/8,  viscous  effects  are  of  the  same  order  as  the 
centrifugal  acceleration  of  the  basic  state  in  the  determination  of  the  Gortler  number,  and  the 
leading  order  inviscid  result  (3.35)  is  to  be  modified.  We  assume  that  to  leading  order  the 
Gortler  number  now  expands  as 


_  2  BM3!2  4 

G  j  \/n  \3/2  4”  ^  SO- 

K{Xn)(  2xn)3/2 


(5.31) 


Here  the  first  term  is  due  to  the  curvature  of  the  basic  state  and  the  second  term  is  due  to 
viscous  effects  and  is  to  be  determined. 

For  convenience,  we  introduce  a  small  parameter  e  and  an  0(1)  constant  N  by 


edM  ~,  N  d='  M3/2c4, 
a 


so  that  (5.31)  can  be  written  as 


C  —  .  —  4- 

«(z„)(2xn)3/2  e4  e4 ' 


(5.32) 


(5.33) 


To  determine  the  higher  order  correction  terms  to  the  Gortler  number  expansion,  we  shall  first 
fix  the  Gortler  number  as  given  by  (5.33)  and  consider  the  evolution  of  Gortler  vortices  in  the 
neighbourhood  of  the  leading  order  neutral  position  xn  defined  by  (5.31),  aiming  at  finding  the 
second  order  correction  say  ex„  to  the  neutral  position.  As  we  have  remarked  at  the  beginning 
of  the  first  subsection,  replacing  x„  by  x„  -  ex„  in  (5.33)  would  give  the  appropriate  expansion 
of  the  Gortler  number  for  vortices  neutrally  stable  at  x  =  x„. 

The  vortices  under  consideration  vary  on  small  lengthscales  in  both  x  and  77  directions. 
In  the  streamwise  direction,  their  growth  rate  can  be  shown  to  be  0(l/e)  so  that  they  evolve 
on  an  0(e)  lengthscale.  In  the  77  direction,  they  are  confined  to  an  0(e1^2)  thin  viscous  layer 
because  of  their  small  wavelength  character.  We  therefore  define  two  new  variables  X  and  £ 

*  =  C  =  (5-34) 


where  77*  is  the  centre  of  vortex  activity  and  is  to  be  determined. 

We  now  look  for  asymptotic  solutions  of  the  form 

t  =  60{X,  0  +  e1/afl1(x,  0  +  MX,  0  +  •  •  • , 

v  =  e~2[Vo(X,C )  +  e1/2Vi(X,()  +  eV2(X,()  +  •••], 
w  =  e-3/2[Wo(X,C)  +  e1/2Wi(X,  Q  +  eW2{X,Q  +  •••], 
P  =  e-5'2[Po(*,0  +  e1!2  Pi(X,Q  +  eP2(X,  ()  +  •••]. 


(5.35) 
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Here  the  relative  orders  of  the  perturbation  quantities  are  deduced  from  the  perturbation  equa¬ 
tions  (3.18)— (3.22).  On  inserting  these  expansions  into  the  perturbation  equations  (3.18)-(3.22), 
expanding  all  coefficients  there  about  x  —  xn  and  77  =  tj*  ,  and  then  equating  the  coefficients 
of  like  powers  of  e,  we  obtain  a  hierarchy  of  equations.  To  leading  order,  the  Gortler  number 
go  in  (5.31)  is  determined  as  a  solvability  condition  for  (Vbj  ^0)  and  is  given  by 


90  ctk0Ti  ’ 


whilst  6q  ,  Wq  and  Pq  are  related  to  V0  by 


/,  _  0T1  rr  _  1  9 V0  n  _  Ho®V0 

2foT*x /°’  V2^TodC  Po  To  d(  ’ 


(5.36) 


(5.37) 


where  T0  =  T( tj*),  Tx  =  fo  =  fi(fQ )  k0  =  /c(xn).  Note  that  (5.36)  is  of  the  same 

form  as  (5.30),  as  we  would  expect. 

To  next  order,  we  obtain  three  expressions  similar  to  (5.37)  for  9X,  Wi  and  Pi  in  terms  of 
V\  and  Vo  and  the  condition  that 

!,=„•  =  0,  (5.38) 


which  implies  that  rj*  is  where  go  attains  its  minimum. 

If  we  carry  on  one  order  higher,  we  find  from  a  solvability  condition  for  (V2,02)  that  Vb 
must  satisfy  the  evolution  equation 


where 


^  _  2(l  +  !)%„3^  _  .(J„o  +  ixVo  =  „ 

a ofj,0  dX 


-  ^0  xn  d*go  1  _ 

=  2xnr02  f  2 BN  (*1  3  Ki  1  \ 

3  IffoKO^Xn)3/2  K0  2ln  K°  2l„  J  ’ 


(5.39) 


(5.40) 


3  l5oKo(2Xn)3'2  «0  2xn  Ko  2xnJ 

and  where  kx  =  «/(xn).  The  solutions  of  (5.39)  which  satisfy  the  conditions  Vo  ->  0  as  |£|  -*  00 
can  be  written  as 


Vo  =  Vom(X,0  =  exp 


3  fob 

4(1  +  <7  )T0xn 


Y  (  . 

X-tm+j)— 


p(-m-i,  (4a)'/4c),  m  =  0,1,2,  (5.41) 

where  U  is  a  parabolic  cylinder  function.  The  neutral  position  x„  can  be  taken  to  be  the  point 
where  dVo/dX  =  0,  so  that  the  mth  mode  is  neutrally  stable  at 


y  -  f  . 

X  =  xn  =  (m  +  j)-j-. 


(5.42) 
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The  most  unstable  mode  corresponds  to  m  =  0.  We  therefore  have 


.  Va 

Xn  ~  T 

With  the  expression  for  xn  determined,  we  can  now  replace  xn  by  in  —  exn  in  (5.33)  and  then 
expand  the  two  terms  on  the  right  hand  side  up  to  and  including  the  0(l/e3)  term,  hence 
obtaining  the  expansion 


2 BM3'2  go  1  1  1 3go  d*g0 

«(*»)(  2xn)3/2  +  e3  '  2f2  Qrjs  +  ’ 


(5.44) 


for  Gortler  vortices  which  are  neutrally  stable  at  position  xn. 

As  is  to  be  expected,  the  Gortler  number  expansion  (5.44)  agrees  in  the  special  case 
k  =  (2x)-3/2  with  the  combination  of  (4.1),  (4.10a),  (4.12d)  and  (4.17)  (note  that  the  small 
parameter  e  there  corresponds  to  e/v/2xn  here).  This  means  that  in  the  large  wavenumber 
limit,  the  relation  (5.44)  is  a  universal  expression  for  the  neutral  Gortler  number,  which  is 
valid  for  all  wall  curvatures. 

Finally  in  this  subsection  we  stress  that  a  more  unstable  version  of  the  parallel  viscous  mode 
can  be  obtained  by  taking  go  bigger  than  the  value  given  by  (5.36),  in  that  case  we  must  allow 
for  a  growth  rate  of  order  M  J  and  then  (5.36)  is  replaced  by  an  equation  to  determine  that 
growth  rate.  This  structure  then  enables  a  direct  connection  between  the  strongly  unstable 
inviscid  mode  and  the  parallel  inviscid  mode  at  relatively  high  Gortler  numbers  with  G-Gn  ~ 
Aft.  The  analysis  for  this  more  unstable  regime  is  essentially  identical  to  that  given  by  Denier, 
Hall,  and  Seddougui  (1990)  in  the  incompressible  case. 


6  The  wall  mode 

It  has  been  established  in  section  3  that  as  the  wavenumber  becomes  large,  Gortler  vortices 
become  increasingly  trapped  in  the  0(1)  temperature  adjustment  layer.  Thus  the  preceding 
three  sections  are  devoted  to  Gortler  vortices  which  have  wavelength  of  0(1)  or  smaller  and 
which  are  trapped  in  the  temperature  adjustment  layer.  Clearly  it  is  possible  for  vortices 
of  wavelength  smaller  than  the  thickness  of  the  transition  layer  to  be  excited,  far  enough 
downstream  the  local  wavenumber  will  become  comparable  to  the  adjustment  layer  thickness 
and  the  previous  analysis  will  apply.  However  before  this  occurs  the  vortices  must  be  described 
by  an  analysis  which  takes  account  of  the  fact  that  they  are  of  wavelength  much  larger  than 
the  adjustment  layer  thickness,  we  shall  now  address  that  situation.  In  fact  it  can  be  deduced 
from  the  definition  (2.11)  and  (3.1c)  that  the  variation  dy  of  the  physical  variable  y  and  the 
variation  of  the  similarity  variable  drj  satisfy 


dy  =  sFIxTdt]. 


(6.1) 


The  wall  layer  which  corresponds  to  r\  =  0{M~ 1/2)  with  t  =  0(M2)  is  therefore  actually  of 
0(M 3/2)  thickness  in  terms  of  the  physical  variable  y,  whilst  the  temperature  adjustment  layer 
is  still  of  0(1)  thickness.  Thus  a  natural  scale  for  larger  wavelength  vortices  is  provided  by  the 
thickness  of  the  wall  layer,  the  appropriate  size  of  the  Gortler  number  is  found  by  rescaling 
the  vortex  wavelength  and  velocity  field  by  the  scales  relevant  to  the  wall  layer.  Such  Gortler 
vortices  are  referred  to  as  the  wall  mode  and  are  studied  in  the  present  section  to  complete 
our  stability  analysis. 

Since  in  the  large  Mach  number  limit  the  boundary  layer  thickens  by  0(M 3/2),  we  should 
rescale  ( y,z )  by  a  factor  M 3/2  and  the  corresponding  velocity  components  likewise.  This 
effectively  replaces  all  “Jl-1/2”s  by  “ii-1/2M3/2”.  It  is  therefore  appropriate  to  rescale  the 
Gortler  number  G  and  the  wavenumber  a  by  defining 

G  =  M~3^2k(x)(2x)3^2G,  k  =  y/2xM3/2a.  (6.2) 

In  the  wall  layer,  the  basic  state  is  from  (2.12)  and  (3.5)  given  by 

AT3/2  - 

u  =  F'(Y),  v  =  - —[-TF(Y )  +  F'(y)fi(r)],  (6.3) 

v  2x 

where 

fi(y)d=  [  TdY.  (6.4) 

Jo 

The  various  partial  derivatives  of  u,v  and  fi  which  appear  in  the  perturbation  equations 
(3.18)-(3.22)  have  to  be  computed  before  we  can  deduce  the  relative  orders  of  the  pertur¬ 
bation  quantities.  Such  expressions  are  given  in  the  Appendix  to  this  paper.  With  the  aid  of 
these  expressions,  we  can  show  from  (3.18)-(3.22)  that  the  relative  scalings  of  the  velocity, 
pressure  and  temperature  disturbance  fields  are  given  by 

V  =  0(M3/2U ),  W  =  0(M3/2£f),  T  =  0{M2U),  P  =  0(MU).  (6.5) 

We  therefore  look  for  solutions  of  the  form 


u  =  ff(*,y)  +  ...,  v  =  M3/2v(x,Y)  +  ---, 


W  =  M3/2W(x,Y)+  •••,  P  =  M2P(x,Y  )  +  •••,  (6.6) 

T  =  M20(x,y)  +  .... 

On  substituting  (6.6)  into  the  perturbation  equations  (3.18)-(3.22)  and  then  equating  the 
coefficients  of  like  powers  of  Af,  we  obtain  to  leading  order  the  following  set  of  partial  differential 
equations  which  govern  the  evolution  of  Gortler  vortices  in  the  wall  layer: 
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(6.11) 
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Obviously,  this  set  of  partial  differential  equations  have  to  be  solved  numerically  to  deter¬ 
mine  the  evolution  properties  of  Gortler  vortices  in  the  wall  layer.  We  further  note  that  the 
downstream  velocity  component  of  the  perturbation  now  does  not  decouple  from  the  other 
disturbance  quantities.  For  the  purpose  of  numerical  calculation,  it  is  convenient  to  eliminate 
the  pressure  perturbation  P  and  W  from  the  above  five  equations.  After  some  manipulation, 
we  obtain 
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The  above  equations  are  to  be  solved  subject  to  the  following  boundary  conditions: 

At  Y  =  0, 

dV 

tf  =  F  =  ^  =  0, 
dd 

=  0,  if  the  wall  is  thermally  insulated  (6.15) 


As  Y  — ►  oo. 


0  =  0,  if  the  wall  is  under  cooling. 


tJ  ->  0,  V  -  0,  ^  -  0,  0  -  0. 


(6.16) 


The  precise  large  Y  decay  behaviour  of  the  perturbation  quantities  can  be  deduced  from 
the  perturbation  equations  (6.12)-(6.14),  but  the  derivation  is  tedious.  However,  the  far 
downstream  limit  of  such  decay  behaviour  can  be  deduced  very  easily  from  the  condition  that 
the  solutions  of  the  wall  mode  under  consideration  in  the  double  limit  k  — >  oo  and  Y  — ►  oo 
should  match  with  the  inviscid  solutions  (3.32)  and  (3.33)  when  k  0  and  77  — ►  0  there.  In  the 
limit  k  — ►  0,  equation  (3.32)  reduces  to  equation  (3.37)  and  the  latter  has  solutions  confined 
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to  a  thin  layer  77  =  0(k 1//3)  near  the  wall.  It  can  be  shown  that  in  the  further  limit  77  —►  0,  the 
solution  of  (3.37)  has 

14-4.  (6.17) 

V 

and  from  (3.33), 

To  -  n't,  -  (6.18) 

where  Ci,ca  and  C3  are  all  functions  of  x.  By  the  above-mentioned  matching  condition,  the 
relations  (6.17)  and  (6.18)  with  77  replaced  by  Y  are  also  the  far  downstream,  large  Y  decay 
behaviour  of  the  perturbation  quantities  in  the  present  problem. 

For  any  initial  disturbance  located  at  a  given  downstream  position,  we  can  determine  the 
position  at  which  the  disturbance  becomes  neutrally  stable  by  integrating  these  equations  using 
the  same  marching  procedure  as  the  one  used  in  solving  the  system  of  equations  (4.3)-(4.7). 
As  an  illustrative  example,  we  impose  the  following  initial  disturbance  at  the  location  x  =  50: 

£f(50,Y)  =  yse-y\  7(50,  Y)  =  0(50, Y)  =  0.  (6.19) 


The  wall  curvature  is  taken  to  be  \/2x  and  the  wall  is  assumed  to  be  thermally  insulated.  For 
a  given  wavenumber  and  a  given  Gortler  number,  we  can  march  downstream  until  we  reach 
the  neutral  position  where  a  certain  energy  measure  has  zero  growth  rate.  We  then  calculate 
the  local  wavenumber  k  and  the  local  Gortler  number  G  at  the  neutral  position.  By  fixing 
the  wavenumber  a(=  k{2x)~x!2)  at  1  and  varying  the  Gortler  number  M~3/2G  from  0.0001  to 
0.03,  we  obtain  a  series  of  neutral  points  ( k,G ).  Fig.13  shows  the  neutral  curves  corresponding 
to  the  following  three  energy  measures: 


Ei  =  (°°  tj2dy  =  V2 xM3/3  [°°  U2TdY, 

Jo  Jo 
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(6.20) 
(6.21) 

jo  -dy'"'-  V2i  Jo  fKdY '  ““  (6,22) 

Fig.14-Fig.17  shows  the  downstream  evolution  of  Gortler  vortices  corresponding  to  the  above 
conditions  with  M~3^2G  =  0.001.  We  observe  that  all  neutral  curves  decrease  monotonically 
with  respect  to  the  local  wavenumber  and  that  Gortler  vortices  become  increasingly  more 
and  more  shifted  to  the  right  (i.e.  towards  the  temperature  adjustment  layer)  as  they  evolve 
downstream.  This  is  certainly  to  be  expected  since  in  the  large  local  wavenumber  limit  the 
wall  mode  has  to  match  with  the  mode  trapped  in  the  temperature  adjustment  layer.  By 
taking  the  large  wavenumber  limit  of  the  wall  mode  equations  we  can  show  that  the  neutral 
curve  should  tend  to  the  limit  6  =  2 B  =  1.2543.  In  order  to  realize  this  limit  numerically  we 
have  to  carry  out  our  calculation  at  very  large  wavenumbers.  This  prc'?nts  some  numerical 
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difficulties,  because  on  one  hand,  if  we  fix  the  Gortler  number  and  vary  the  wavenumber, 
the  initial  conditions  soon  become  incompatible  with  the  differential  equations  since  for  large 
wavenumbers  Gortler  vortices  have  to  be  trapped  in  a  region  away  from  the  wall;  on  the 
other  hand,  if  we  fix  the  wavenumber  and  vary  the  Gortler  number,  a  large  local  wavenumber 
corresponds  to  a  large  downstream  neutral  position  and  there  contamination  from  the  “finite 
infinity”  becomes  important  because  of  the  algebraic  decay  behaviour  (implying  that  we  have 
to  choose  a  larger  infinity).  Finally,  we  remark  that  the  2-derivative  of  6  in  the  expressions  for 
W3  and  W  could  be  eliminated  with  the  aid  of  equation  (6.13)  to  give  a  different  formulation 
and  thus  to  provide  a  check  on  our  numerical  scheme.  We  have  done  so  and  have  obtained 
identical  results. 

We  conclude  this  section  by  stating  the  most  important  results  of  our  investigation  of  the 
wall  mode.  We  have  shown  that  the  wall  layer  can  support  a  disturbance  trapped  in  the  wall 
layer  with  wavelength  comparable  with  the  wall  layer  thickness.  This  mode  is  dominated  by 
nonparallel  effects  and  has  a  neutral  Gortler  number  which  is  a  monotonic  decreasing  function 
of  the  vortex  wavenumber.  In  the  limit  of  high  vortex  wavenumbers  the  mode  takes  on  a  struc¬ 
ture  essentially  identical  to  that  found  for  the  small  wavenumber  limit  of  the  inviscid  modes  of 
wavelength  comparable  with  the  adjustment  layer  thickness.  Moreover  in  this  limit  the  vortex 
has  a  neutral  Gortler  number  which  approaches  from  above  the  zeroth  order  approximation  to 
the  neutral  Gortler  number  of  the  wall  modes. 

7  Real  gas  effects 

In  our  previous  discussions,  we  have  assumed  that  the  fluid  under  consideration  is  an  ideal  gas 
undergoing  no  dissociation.  Our  asymptotic  analysis  based  on  the  large  Mach  number  limit 
has  yielded  some  revealing  results  about  the  stability  properties  of  hypersonic  boundary  layers. 
However,  in  the  large  Mach  number  limit,  we  would  expect  that  the  wall  temperature  should 
be  well  above  the  temperature  at  which  dissociation  takes  place.  Take  a  boundary  layer  over 
a  thermally  insulated  wall  as  an  example.  The  temperature  at  the  wall  is  given  by 

Tw  =  i(7  -  \)M2T„ 

when  the  Prandtl  number  is  unity.  At  a  standard  altitude  of  53  km,  the  air  temperature  T «, 
is  283K.  If  we  take  7  =  1.4,  M  =  25,  then  Tw  =  35,375  K.  Since  at  the  far  lower  temperature 
of  2500  K,  the  oxygen  molecules  in  the  air  have  already  begun  to  dissociate,  it  is  clear  that  an 
investigation  which  takes  gas  dissociation  into  account  is  vitally  important! 

A  complete  theory  on  real  gas  effects  should  at  least  incorporate  the  following  important 
chemical  reactions: 

02  -  20,  N2  **  2 N, 
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N  +  O^NO,  N  +  0  ^  N0+  +  e~. 


However,  in  order  to  expose  the  most  important  features  and  at  the  same  time  keep  the  algebra 
to  a  minimum,  we  begin  our  investigation  by  eliminating  the  less  essential  complications  which 
arise  from  the  detailed  composition  of  air  and  confining  our  attention  to  a  pure  dissociating 
diatomic  gas.  The  dissociation  process  is  therefore  denoted  by 


A  +  A. 


(7.1) 


In  our  following  discussion,  02  and  1V2  will  be  chosen  for  numerical  illustration.  F" 
we  assume  that  the  fluid  under  consideration  is  an  ideal  dissociating  gas  which 
relation 
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(see  Becker  (1968),  p.  36,  or  Lighthill  (1957),  p.6).  Here  a,  p  and  T  have  the  same  meanings 
as  defined  in  §2,  whilst  pd  and  Td  are  respectively  the  characteristic  pressure  and  temperature 
for  dissociation. 

For  the  gas  mixture  of  A  and  A 2,  Dalton’s  law  gives 


p  =  pi 


7 ti23?T 
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(7.3) 


V  V 

where  pi  is  the  pressure  which  the  ith  component  would  exert  individually  if  alone  in  volume 
V  at  temperature  T,  ni  are  mole  numbers  and  Si  is  the  universal  gas  constant.  Here  we  use 
subscript  1  to  signify  component  gas  A2  and  subscript  2  for  A.  Assume  that  the  weight  of  a 
mole  of  gas  A  is  m,  the  molar  weight  of  gas  A2  is  then  2m.  We  therefore  have 
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With  the  use  of  (7.4)  and  (7.5),  we  can  rewrite  the  equation  of  state  (7.3)  as 


(7.4) 

(7.5) 


p  =  (1  +  a  )$lpT, 


(7.6) 


where  &  =  Si/(2m)  is  a  gas  constant. 

On  inserting  (7.6)  into  (7.2),  we  obtain 

J?_  =  Sie-TtIT  (7.7) 

1  -  a  p 

where  pd  =  Pd/(^tTd)  is  the  characteristic  density  for  dissociation.  pd,  Pd  and  Td  are  in 
general  functions  of  T,  but  their  variations  over  a  large  temperature  range  are  very  slight. 
Their  typical  values  are  given  in  Table  2  (taken  from  Lighthill  (1957)). 
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Table  2 


Td(°K) 
Oxygen  59000 

Nitrogen  113000 


Pd(atoms)  pd(gm/cm3) 

2.3  x  107  150 

4.1  x  107  130 


It  should  be  noted  that  for  atmospheric  values  of  p,  pd/ p  is  at  least  105.  Thus  although 
pd,  Td  and  pd  are  called  characteristic  quantities,  they  are  not  typical  of  the  actual  values 
of  p,  T  and  p  at  all.  For  pdj p  =  105,  (7.7)  shows  that  a  is  0.05  (5%  dissociation)  when 
T /Td  =  0.057,  and  a  is  0.95  (95%  dissociation)  when  TJTd  =  0.116.  For  densities  typical  of 
the  upper  atmosphere,  with  (say)  pd/ p  =  107,  these  values  of  T jTd  would  be  reduced  to  0.045 
and  0.076  respectively.  For  a  fixed  value  of  p^/p  =  4  x  108  and  Td/Too  =  4  x  103,  the  variation 
of  a  with  respect  to  T/T^  is  plotted  in  Fig. 18. 

Relation  (7.7)  was  first  obtained  by  Lighthill  (1957)  from  quantum  mechanics,  and  there¬ 
fore  the  ideal  dissociating  gas  discussed  here  is  also  called  the  Lighthill  gas.  Using  a  purely 
mathematical  argument,  Becker  (1968)  has  shown  that  the  Lighthill  gas  is  a  special  case  of  a 
more  general  class  of  gases.  An  important  property  of  Lighthill  gases  is  that  A2  and  A  have 
the  same  specific  heats  at  constant  volume,  that  is, 


Cy  1  =  Cv2- 


(7.8) 


7.1  Constitutive  properties  of  a  dissociating  gas 

In  order  to  see  the  complications  which  arise  from  gas  dissociation,  let  us  first  recall  that  an 
ideal  gas  has  the  following  properties  which  have  greatly  simplified  our  previous  analysis: 

(i)  The  equation  of  state  has  the  simple  form  P  =  ftpT  as  compared  with  (7.6); 

(ii)  Specific  heats  cp  and  c„  can  be  taken  to  be  constant; 

(iii)  The  coefficient  of  thermal  conduction  k  is  proportional  to  the  shear  viscosity  p  and  the 
Prandtl  number  a  =  pcp/k  is  usually  taken  to  be  constant; 

(iv)  The  viscosity  can  be  taken  to  be  related  to  the  temperature  by  Chapman’s  law,  or  more 
accurately,  by  Sutherland’s  law. 

When  part  of  an  ideal  gas  has  been  dissociated,  the  gas  becomes  a  mixture  of  two  component 
gases  and  all  of  the  above  properties  are  changed.  We  have  shown  above  that  the  equation  of 
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state  is  modified  to  the  form  given  by  (7.6).  In  this  subsection,  we  derive  the  corresponding 
expressions  for  (x,  k  and  cp.  Our  derivation  is  based  on  work  discussed  in  “ The  mathematical 
theory  of  non-uniform  gases ”  by  Chapman  and  Cowling  (1970). 


The  viscosity  n 

In  the  kinetic  theory  of  gases,  different  expressions  for  the  transport  properties  fx  and  k  have 
been  found  depending  on  what  model  is  used  for  the  interaction  of  the  gas  molecules.  For 
example,  for  a  simple  gas  which  consists  of  smooth  rigid  elastic  spherical  molecules,  fx  is 
proportioned  to  VT,  whilst  if  the  gas  is  taken  to  consist  of  smooth  rigid  elastic  spherical 
molecules  each  of  which  is  surrounded  by  a  field  of  attractive  force,  fx  is  given 

5 


(7.9) 


where  c  is  the  diameter  of  the  molecule,  5  the  potential  energy  of  the  mutual  attraction  of  two 
molecules  when  in  contact,  m  the  molecular  weight.  Relation  (7.9)  is  known  as  Sutherland’s 
viscosity  law. 

For  a  binary  mixture  of  gases,  Wilke’s  law  gives  an  approximate  expression  for  the  viscosity 
of  the  mixture: 

(7.10) 


fx  = 


Xi 


X3 


aii a. 
Mi  Mia 


_Ei_  _l  ’ 
Mai  +  Ma 


where  ^  and  /z2  are  the  viscosities  of  the  two  component  gases;  xx  and  x2  denote  the  propor¬ 
tions  by  volume  of  the  two  gases  in  the  mixture.  Since  by  Avogadro’s  law  equal  volumes  of 
different  gases  contain  an  equal  number  of  molecules  or  moles,  we  have 


7li  i  -  u  1X2  ,7„< 

xj  = - -  - - ,  x2  = -  =  — — .  (7.11) 

71 1  +  712  1  +  Ct  Tl\  +  712  1  4"  OL 

In  (7.10),  /xj2  and  fx2i  are  the  mutual  viscosities  of  the  two  component  gases  and  are  given  by 

[  KT(m!  +m2)l1/2 


Q 


n2 


2  a 


M12  - 


5mi 

16cf2 


/(1  +  ^r),  M21  =  — -M12  -  xMn.  (7-12) 

1  77li  Z 


2rmi  771 2 

where  c12  can  be  taken  to  be  (ci  +  c2)/ 2,  whilst  the  expression  for  5x2  has  to  be  found  empir¬ 
ically.  Lindsay  and  Bromley  (1950)  suggested  that  5i2  =  \/^i^2-  Here  and  Sx  (i=  1 ,2)  are 
constants  appearing  in  (7.9),  associated  with  fX{. 

With  the  add  of  (7.9)  and  (7.12),  we  can  now  write  down  the  appropriate  expressions  for 
the  viscosities  of  the  component  gases.  These  are 


fXi  = 


Here 


A,r3/2 

T  +  5,’ 


= 


i  =  1,2,  /x12  = 


a3t3/2 


I  SRm, 


<43  = 


16c?  V  * 

where  we  have  made  use  of  the  fact  that  mi  =  2m2 


M21  =  M12/2, 


(7.13) 


(7.14) 
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Coefficient  of  thermal  conduction  k 

For  a  pure  gas,  the  coefficient  of  thermal  conduction  is  related  to  the  viscosity  p  by  the  simple 
formula 


k  —  fcvfi, 

where  /  is  a  constant  and  7//  is  usually  defined  as  the  Prandtl  number. 
For  a  binary  gas  mixture,  Wassiljewa’s  formula  gives 

k  =  Xl  +  Xi 


(7.15) 


(7.16) 


where  k\  and  k 2  are  the  thermal  conductivities  of  the  two  component  gases  and  Aj12,  and  k2i 
are  the  mutual  thermal  conductivities. 


It  is  suggested  in  Chapman  and  Cowling  (1970,  p.256)  to  put 


k\  _  Du 
ku~ai2Di2' 


k2  D22 

7—  =  7T-  > 

«21  D12 


(7.17) 


where  the  coefficients  07 2  and  a2i  are  regarded  as  functions  of  mijm2,  which  are  determined 
semi-empirically,  and  for  Sutherland’s  model  the  diffusion  coefficients  Du,  D 12  and  D22  are 
given  by 

6  6  6 

Du  =  Z - Mil  -®12  =  r - Pl2i  D22  =  T - P-2,  (7.18) 

bpi  5pi  bp2 


where  pi  is  the  density  of  the  ith  component  gas,  when  pure,  at  the  pressure  and  temperature 
of  the  actual  mixture. 

If  we  take  ai2  =  a2i  =  1  for  simplicity  and  further  make  use  of  the  above  relations,  relation 
(7.16)  then  becomes 

l- .  ,  \  Jl  .  (A/A)<*»  1 


-c„/i{l 


-l  1  1  xla.ii 

n  Mia  ia  Mia 


(7.19) 


where  we  have  used  the  fact  that  for  an  ideal  dissociating  gas,  c„i  =  c„2  —  c„;  whilst  f\  and 
f2  are  the  constants  appearing  in  (7.15)  corresponding  to  gases  A2  and  A  respectively. 

For  easy  comparison,  we  rewrite  (7.10)  here,  making  use  of  the  relation  (7.12), 


14 .  2Z  .  i£L  ~r  l4.2iL.ii2_- 
xi  Mia  '  *a  Mia 


(7.20) 


Comparing  (7.19)  with  (7.20)  shows  that  we  can  almost  write  k  =  cvfip  for  the  gas  mixture, 
only  if  f2  —  fi-  It  is  therefore  desirable  to  investigate  the  values  of  fi  (for  diatomic  gases)  and 
/2  (for  monatomic  gases). 

For  all  smooth  spherically  symmetrical  molecules,  it  has  been  shown  that  taking  f2  =  2.5 
provides  a  very  good  approximation.  However,  for  diatomic  and  polyatomic  gases,  a  variety  of 
expressions  have  been  suggested  for  /.  One  of  these  is  Eucken’s  formula: 


f  =  \{  97-5). 


(7.21) 
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This  formula  is  valid  under  the  assumption  that  the  transport  of  momentum  and  translational 
energy  is  unaffected  by  the  internal  molecular  motions,  and  that  internal  energy  is  transported 
at  the  same  rate  as  momentum. 

For  diatomic  gases,  there  are  3  translational  degrees  of  freedom,  two  rotational  degrees  of 
freedom  and  on  the  average  one  vibrational  degree  of  freedom  (may  vary  between  zero  and 
two  depending  on  whether  vibration  has  been  fully  excited  or  not).  Thus  the  internal  energy 
is  given  by 


1 

€l  ~  2 


3?T  .  n  ,,  33 ?T 

P^)(3  +  2+1)-'2^T> 


(Recalling  that  2 m  is  the  molecular  weight  of  *4).  Consequently,  the  specific  heats  are 


dt\  33? 


3 f? 


We  then  have 


Cvl  “  BT  -  2 m*  ^  ~  Cvl  +  2m 


^  =  ^H  =  133’  f=l  =  lJ5- 


(7.22) 


(7.23) 


(7.24) 


At  lower  temperatures  when  the  vibrational  mode  is  not  excited, 

5.8?T,  5.8?.  7.3?. 

61  ~  2  2m  Cul  —  2  2m  1  ^  -  2^2m^’ 


and  therefore 


/  =  1.90. 

At  the  other  limit  when  the  vibrational  mode  is  fully  excited, 


7  =  -  =  i.4, 


and  hence 


7.3?T,  7.3?,  9.3?, 

Cl  ~  2  2m  ’  Cvl  _  2  2m  Cpl  _  2^2m^’ 


7=7  =  1-29,  /  =  1-68- 


Therefore,  in  general, 

1.68  <  /  <  1.90.  (7.25) 

In  passing,  we  note  that  for  monatomic  gases,  there  are  only  three  translational  degrees  of 
freedom  and  therefore 

1  q  3Rr 

e2  =  -  •  —  ■  3  = -T— .  (7.26) 


Comparing  (7.26)  with  (7.22)  then  shows  that  monatomic  and  diatomic  gases  have  internal 
energies  differing  only  on  their  zero  point  energies  and  that  the  basic  assumption  cvi  =  cv 2  for 
an  ideal  dissociating  gas  is  indeed  valid.  Also,  the  /  value  2.5  for  a  monatomic  gas  can  actually 
be  read  off  from  Eucken’s  formula  (7.21)  by  noting  that  cv 2  =  38?/(2m),  cP2  =  5S?/(2m)  and 
7  =  5/3. 
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We  now  return  to  (7.19).  From  the  above  discussion,  it  is  appropriate  to  put  cv  = 
3R/(2m),  /i  =  7/4,  f2  =  5/2.  The  coefficient  of  thermal  conduction  k  for  the  mixture 
is  then  given  by  the  formula 


_  21  ft  f  Mi  .10  M2  1 

8m  11  +  7  '  1  +  J  ' 

v.  Mis  1  ij  Mia  ' 

Specific  heats  of  the  mixture 

The  specific  internal  energy  for  A  and  A2  both  have  the  same  expression 


(7.27) 


The  total  internal  energy  of  the  mixture  is  then  given  by 

E  =:  — — (2m  •  Tii  +  m  ■  712). 
2m 


(7.28) 


Since  A  and  A2  do  not  necessarily  have  the  same  zero-point  energy,  relation  (7.28)  should  be 
modified  to 

1 

E  —  — — (2m  •  ni -)- m  •  712)  + -3Wd«2,  (7.29) 

2m  2 

where  the  last  term  represents  the  difference  of  zero-point  energies  of  A  over  A2- 

On  dividing  (7.29)  by  the  total  mass  of  the  mixture  (2mni  +mn2),  we  obtain  the  expression 
for  the  specific  internal  energy  of  the  mixture: 


e  =  3&T  +  m’da. 


(7.30) 


With  the  aid  of  this  relation  and  the  equation  of  state  (7.6),  we  can  easily  calculate  the 
specific  enthalpy  as  follows: 


h  =  e  +  Z  =  (  4  +  a)&T  +  %Tda. 


(7.31) 


To  calculate  specific  heats,  we  have  to  first  of  all  evaluate  ( da/dT)p  and  ( da/dT)p .  From 

(7-7), 

<§?>'  =  <&v  =  £(!  +  £)■ ■  s-f1  ~ “*>•  (732) 

It  then  follows  from  (7.30)  and  (7.31)  that 

c.  =  (^).=3S  +  S(|!)!.2^f^,  ct=(4  +  a)S+|(l  +  ^)1.a(l-a:).  (7.33) 

The  quantity  ( dh/dp)T  which  appears  in  (2.8)  can  now  be  calculated  with  the  use  of  (7.31) 
and  (7.2).  The  result  is 


,on  $t(T  +  Td)  2\  J-  ,  x 

(^)t  = - Yv —  ’  a(l  -  °  )  =  -~(1  +  tt)  •  a\l  -  °0- 


— (1  +  — ' 
2 pK  r  T‘ 


(7.34) 
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7.2  Modification  of  the  basic  state 

As  we  have  remarked  before,  the  basic  state  equations  (2.13)  and  (2.14)  are  independent  of 
constitutive  assumptions  and  therefore  also  valid  for  the  gas  mixture  under  consideration, 
although  now  p,  p,  k  and  Cp  are  calculated  from  the  more  complicated  expressions  (7.6), 
(7.10),  (7.27)  and  (7.33),  respectively.  In  these  expressions,  the  function  a  is  given  by  (7.7). 
After  non-dimensionalization,  these  equations  become 


where 


9  (l  +  a)r 

(1  +  m1)f3/2  izfl  +  mQf3/2 

f  +  +  i  '  T^(f  +  ™3)  T  +  rhi  +  ^”(r  +  m3)’ 

(l  +  mi)f3/2  |  (10ia/7)(l  +  mQf3/2 

f  +  rhi  +  ^  •  +  m3)  f  +  rh2  +  •  l=a(f  +  m3) ' 

^ = 1 + + s(i + o' ■  “(i  _ 

1  -  o2  Poo 


5i  52  _  5i5a 

mi  =  — ,m2  =  =-,m3  =  -7= — , 

■*  00  ■*  00  00 


(7.35) 


(7.36) 


(7.37) 


(7.38) 


(7.39) 


A2  =  4i  = 


7  _  A3  _  /3  c? 

;  3  “  Ai  V  2  c2,' 


(7.40) 


J  Ax  V2c2’  °  Ax  V  2  c2a  ' 

For  the  purpose  of  asymptotic  analysis,  it  is  convenient  to  define  two  new  constants  a  and  b 

°  =  (741) 


Then  equation  (7.39)  becomes 


or  aT 

_  =  _ 

1  -  a2  M2 


(7.42) 


which  displays  the  physical  fact  that  dissociation  will  take  place  in  the  hottest  region  where 
f  =  0(M2). 

An  asymptotic  analysis  of  equations  (2.13)  and  (2.14)  shows  that  the  boundary  layer  struc¬ 
ture  in  the  present  case  is  similar  to  that  for  an  ideal  gas.  In  particular,  the  boundary  layer 
can  be  divided  into  two  regions:  an  inner  region  77  =  0(M-1/2)  and  a  temperature  adjustment 
region  77  =  0(1).  In  the  inner  region,  we  define  new  variables  V,  T  as  in  (3.5),  but  now  (3.6) 
and  (3.7)  are  replaced  by 


(7.43) 
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where 


1  +  mg 
a 


t'  y 

Vf) 


+  CpFT'  +  (7  -  1)(1  +  mg) 


hi(a)  {F"f 
1  +  a  s/f 


=  0, 


(7.44) 


a,  M  -  °)  + 


A2A3a 


2a  -f  A3(l  -  ol)  A3a  +  A2(  1  -  a)  ’ 


•  aci 

flo  ~ 


def.  A3(l  -  a) 


+ 


(10A2A3/7)a 


2a  +  A3(l-a)  A3a  +  A2(l-a)’ 


cp  =  l  +  ia  +  i(l  +  |)2.a(l-a2)) 


(7.45) 

(7.46) 

(7.47) 


and  where  a  satisfies 

JZT~ 2  =  ^e~blf.  (7.48) 

As  Y  -»  00,  T  -*  0  and  we  expect  that  a  -*  0.  Prom  (7.45),  (7.46)  and  (7.38)  we  have 
hg  ►  1,  h2  *  1  and  £p  — >  1.  Equations  (7.43)  and  (7.44)  can  then  be  approximated  by 


(I  +  tth)^)'  +  FF"  =  0, 
VT 


1  +  mg  T' 

^  W 

which  have  asymptotic  solution 

D 


F  =  P-/3  + 


( Y  -  P)3/* 


+  •••,  f  = 


3(1  +  mg) 


(7.49) 

(n2  n 

Vt  ’ 

(7.50) 

la  1 

(Y-P)4+  •’ 

(7.51) 

where  both  D  and  P  are  constants.  Therefore,  equations  (7.43)  and  (7.44)  are  to  be  solved 
subject  to  the  boundary  conditions  (3.8)  and  the  asymptotic  conditions  (7.51).  Results  from 
such  a  numerical  integration  are  shown  in  Fig.l  together  with  those  results  corresponding  to 
the  undissociated  model. 

In  the  temperature  adjustment  layer  77  =  0(1),  a  is  exponentially  small  which  means  that 
no  dissociation  takes  place  in  this  region.  The  basic  state  equations  (2.13)  and  (2.14)  reduce 
to  (3.13)  and  (3.14)  which  are  appropriate  for  an  ideal  gas  and  which  are  to  be  solved  subject 
to  the  same  matching  conditions  (3.15)  and  the  conditions  (3.16)  at  infinity. 


7.3  Modification  of  the  stability  properties 

Since  the  boundary  layer  structure  is  similar  to  that  for  an  ideal  gas  when  dissociation  is  taken 
into  account,  the  qualitative  stability  properties  are  also  similar  with  appropriate  quantitative 
modifications.  In  particular,  we  can  still  show  that  the  mode  trapped  in  the  temperature  ad¬ 
justment  layer  is  most  susceptible  to  Gortler  vortices  and  the  neutral  Gortler  number  expands 
as  in  (3.35),  but  now  the  coefficient  B  is  modified  since  it  is  the  leading  order  contribution  to 
the  integration  of  basic  state  temperature  across  the  whole  boundary  layer.  As  for  higher  order 
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correction  terms  to  the  Gortler  number  expansion,  The  effect  of  gas  dissociation  is  dependent 
on  the  wall  curvature.  In  the  special  curvature  case,  since  the  second  term  G  in  (4.1)  is  de¬ 
termined  by  solving  partial  differential  equations  in  the  undissociated  temperature  adjustment 
layer,  it  is  not  affected  by  consideration  of  gas  dissociation;  whilst  in  the  more  general  curva¬ 
ture  case,  we  can  see  from  (5.9)  that  gas  dissociation  affects  the  second  order  correction  term 
through  E,  but  in  the  higher  wavenumber  case  discussed  in  §5.3,  the  second  order  correction 
term  is  not  affected  by  gas  dissociation,  as  is  clear  from  (5.44). 

In  Fig.19,  we  have  plotted  the  variation  of  B  with  respect  to  the  cooling  coefficient  both  for 
the  real  gas  model  discussed  here  and  for  the  ideal  gas  model  used  in  section  3.  It  is  clear  from 
Fig.19  that  the  values  of  B  is  decreased  by  gas  dissociation  as  well  as  by  wall  cooling.  Therefore, 
both  gas  dissociation  and  wall-cooling  axe  destabilizing.  In  our  numerical  integration  of  the 
boundary  layer  equations  (7.43)  and  (7.44),  we  have  taken  rhi  =  fh=  0.508,  Aj  =  ^3  =  1,  o  = 
1.225  x  109,  6  =  3.2653.  We  have  also  repeated  our  calculation  for  a  few  sets  of  different  values 
for  the  above  five  constants.  We  find  that  the  above  prediction  is  still  valid. 

In  order  to  determine  the  effects  of  gas  dissociation  on  the  wall  mode,  we  have  integrated 
the  perturbation  equations  (6.12)— (6.14)  with  the  basic  state  given  by  the  solutions  of  (7.49) 
and  (7.50),  subject  to  the  same  conditions  as  those  used  to  produce  the  neutral  curves  in 
Fig.13.  Fig.20  to  Fig.22  give  a  comparison  of  the  neutral  curves  corresponding  to  the  ideal  gas 
model  and  the  real  gas  model  with  dissociation.  We  observe  that  in  each  of  these  Figures,  the 
two  neutral  curves  corresponding  to  the  two  models  intersect.  Therefore,  gas  dissociation  can 
have  either  a  destabilizing  effect  or  a  destabilizing  effect  on  the  wall  mode. 

8  Further  discussion 

We  have  seen  above  that  the  Gortler  mechanism  in  a  hypersonic  boundary  layer  of  a  Sutherland 
law  fluid  behaves  in  a  predominantly  viscous  or  inviscid  manner  depending  on  whether  or  not 
the  wall  curvature  varies  like  a: "3  where  x  denotes  distance  along  the  wall.  When  the  wall 
curvature  does  not  have  this  special  form  the  vortices  evolve  over  almost  the  whole  of  the 
wavenumber  space  in  a  nonparallel  manner  subject  to  viscous  effects.  The  only  exception  to 
this  case  is  at  extremely  high  wavenumbers  where  the  vortices  evolve  in  a  quasi-parallel  manner 
essentially  identical  to  that  described  for  incompressible  flows  by  Hall(1982).  Since  the  special 
curvature  distribution  is  possibly  of  little  physical  relevance  we  shall  now  concentrate  on  the 
results  we  have  found  for  the  more  general  curvature  situation. 

The  results  we  have  found  for  the  different  wavenumber  regimes  for  the  general  curvature 
distribution  are  summarized  in  Table  3  below. 
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Table  3 


wavenumber 

wall  mode  0(M-?) 

strongly  unstable  inviscid  mode  0(1) 
near  neutral  inviscid  mode  0(1) 

nonparallel  viscous  mode  O(Afi) 

parallel  viscous  mode  0(M?) 


growth  rate 
0(1) 
O(MT) 

0(mI) 

0(M3) 

O(Jlfi) 


neutral  value  of  G 

O(Aft) H - 

no  neutral  G 
Gn  +  O(M)  +  ••• 

Gjv  +  0(M)  +  ... 
0(M§)  +  ... 


where  Gn  is  as  defined  below  (3.35). 


The  wall  mode  which  was  discussed  in  section  6  is  in  fact  the  counterpart  of  the  so-called 
‘acoustic  mode’  of  the  inviscid  hypersonic  instability  theory  of  flat  plate  boundary  layers.  The 
latter  Rayleigh  instability  has  been  discussed  by  Cowley  and  Hall(1990)  and  Smith  and  Brown 
(1990)  for  Chapman  law  fluids,  and  a  limited  discussion  of  this  mode  for  Sutherland  law  fluids 
can  be  found  in  Blackaby,  Cowley  and  Hall(1990).  The  acoustic  mode  and  the  wall  mode 
discussed  in  section  6  have  the  property  that  they  are  concentrated  in  the  wall  layer  where  the 
streamwise  velocity  component  of  the  basic  state  varies  from  zero  at  the  wall  to  almost  it’s  free- 
stream  value.  The  Rayleigh  acoustic  mode  in  general  evolves  in  a  quasi-parallel  manner  though 
in  the  presence  of  strong  shocks  this  is  not  necessarily  the  case.  The  wall  mode  discussed  in 
this  paper  evolves  in  a  nonparallel  manner  and  becomes  progressively  concentrated  towards 
the  edge  of  the  wall  layer  as  the  local  wavenumber  increases. 

When  the  wavenumber  of  the  vortex  becomes  0(1)  then  the  disturbance  modifies  itself  so 
as  to  become  concentrated  in  the  adjustment  layer  where  the  basic  state  temperature  adjusts 
rapidly  to  it’s  free-stream  value.  The  counterpart  of  this  mode  in  the  Rayleigh  instability 
problem  is  the  so-called  ’vorticity’  mode  investigated  by  Blackaby,  Cowley  and  Hall  (1990)  for 
Sutherland  law  fluids,  and  by  Smith  and  Brown  (1990)  for  Chapman  law  fluids.  We  found  in 
the  present  paper  that  when  the  Gortler  number  G  is  as  given  by  (3.27)  with  G*  >  Q  and  Q 
given  by  (3.28),  then  the  appropriate  expansion  of  the  disturbance  field  is  given  by  (3.31).  Thus 
the  mode  has  growth  rate  0(Afi )  and  the  growth  rate  is  given  by  the  solution  of  the  eigenvalue 
problem  specified  by  (3.32)  subject  to  the  condition  that  Vo  should  vanish  at  0,oo.  Fig.  3 
shows  that  the  growth  rate  increases  monotonically  from  zero  as  the  wavenumber  increases 
and  tends  to  a  constant  at  large  wavenumbers.  At  large  wavenumbers  the  growth  rate  can  be 
matched  onto  the  parallel  viscouB  mode  growth  rate  when  the  vortex  wavenumber  becomes 
formally  O(Aft).  We  note  here  that  the  inviscid  mode  matches  directly  onto  the  parallel 
viscous  mode  structure  at  sufficiently  high  Gortler  numbers.  This  is  exactly  the  situation  with 
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the  inviscid-viscous  connection  between  temporally  growing  Gortler  vortices  in  incompressible 
flows  and  is  a  direct  consequence  of  the  fact  that  the  growth  rate  of  the  inviscid  mode  shown 
in  Fig.3  tends  to  a  constant  at  large  values  of  the  wavenumber.  In  the  corresponding  spatially 
growing  problem  for  incompressible  flows  Denier,  Hall  and  Seddougui  (1990)  show  that  an 
intermediate  region  is  required  to  make  the  inviscid-viscous  connection  and  indeed  that  the 
maximum  possible  growth  rate  is  achieved  in  that  interval.  We  stress  that  this  intermediate 
regime  has  no  connection  with  the  nonparallel  viscous  mode  discussed  in  section  5.  The  latter 
regime  is  appropriate  only  to  the  inviscid-viscous  connection  problem  at  mildly  supercritical 
values  of  G.  In  Fig.  23  we  have  sketched  the  dependence  of  the  growth  rate  on  the  the 
wavenumber  for  the  cases  when  the  Gortler  number  is  mildly  supercritical  (ie  when  G  differs 
only  slightly  from  its  value  required  to  overcome  the  strong  0(M%)  curvature  of  the  basic 
state)  and  the  strongly  unstable  case  with  G  ~  0(M 5 ).  The  broken  parts  of  the  curve  denote 
regimes  where  the  growth  rate  is  given  by  a  nonparallel  calculation.  In  Fig.24  we  sketch  the 
neutral  curve  in  the  local  Gortler  number-local  wavenumber  plane.  In  this  figure  we  have 
indicated  the  corrections  to  the  neutral  curve  associated  with  nonparallel  effects. 

There  are  no  available  experimental  results  with  which  we  can  compare  with  our  calcula¬ 
tions,  this  is  because  experiments  at  the  high  Mach  numbers,  say  10  —  30,  appropriate  to  our 
work  are  exceedingly  difficult  to  perform.  Thus  for  design  purposes  it  is  fair  to  say  that  a 
theoretical  approach  is  the  only  means  at  this  stage  to  predict  the  likely  evolution  of  Gortler 
vortices  in  growing  hypersonic  boundary  layers.  Current  transition  prediction  methods  are  all 
based  on  some  amplitude  growth  criterion  based  on  the  linear  growth  of  a  disturbance.  We 
have  seen  above  that  for  a  realistic  hypersonic  boundary  layer  the  regime  where  vortex  growth 
is  likely  to  occur  is  the  strongly  unstable  inviscid  one.  Hence  in  any  transition  prediction 
method  for  a  hypersonic  boundary  layer  it  would  be  appropriate  to  simply  compute  the  local 
Gortler  number  and  obtain  the  corresponding  growth  rate  from  Fig.3;  this  of  course  should 
only  be  done  for  wavenumbers  less  than  the  neutral  value  associated  with  parallel  viscous  mode 
structure. 

Finally  we  close  by  making  a  few  remarks  about  the  results  of  section  7  which  concerned  the 
effect  of  gas  dissociation  a..«.i  wall  cooling  on  vortex  growth.  First  let  us  make  a  few  remarks 
about  the  effect  of  wall  cooling,  since  the  most  unstable  vortices  correspond  to  the  strongly 
unstable  inviscid  mode  we  concentrate  on  that  situation.  In  that  case  the  main  effect  of  wall 
cooling  or  real  gas  effects  is  to  alter  the  quantity  B  which  fixes  the  scaled  Gortler  number 
above  which  instability  can  occur.  Fig.  19  shows  that  for  an  ideal  gas  the  effect  of  wall  cooling 
is  to  reduce  B  by  a  factor  of  about  .5  when  the  wall  temperature  is  cooled  from  it’s  adiabatic 
value  by  a  factor  of  10.  We  find  that  the  eigenvalue  problem  associated  with  (3.32)  has  the 
maximum  value  of  (3  altered  only  by  a  small  amount  when  this  happens  so  that  at  a  given 
value  of  the  Gortler  number  the  growth  rate  is  increased  by  a  factor  of  about  1.4  when  the 
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wall  temperature  is  reduced  by  a  factor  of  10  from  it’s  adiabatic  value.  By  contrast  real  gas 
effects  have  a  negligible  impact  on  the  growth  rate  of  the  strongly  unstable  mode.  Thus  for 
example  in  Fig.  19  we  see  that,  at  a  fixed  value  of  the  cooling  coefficient,  B  varies  by  only 
about  10  percent  when  the  real  gas  model  is  used.  Thus  the  critical  Gortler  number  or  the 
disturbance  growth  rate  is  altered  only  slightly  by  real  gas  effects  so  it  would  seem  that  any 
transition  prediction  method  could  quite  sensibly  ignore  such  a  complication;  certainly  the 
error  associated  with  doing  so  would  be  negligible  compared  with  the  inherent  error  of  the 
prediction  method. 
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Fig.1  Basic  state  temperature  in  the  wait  layer  corresponding  to 
two  models.  Here  n  Is  the  cooling  coefficient. 
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Fig. 2.  Basic  state  temperature  in  the  adjustment  layer. 
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>/2i(i2l{G'  -  Q) 


Fig.  3  The  dependence  of  the  eigenvalue  y/2x(32/(Gm  -  Q)  of  (3.32)  on  the 
local  wavenumber  k,  which  agrees  with  the  results  from  asymptotic  analysis 
that  as  it  — ♦  0,  y/2x(32/(Gm -Q)  is  of  order  k;  whilst  as  k  -*  oo,  \f2x0i1  j{G*  - 
Q)  tends  to  the  constant  0.5786. 
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Fig.  4  The  eigenfunctions  of  (3.32)  corresponding  to  a  set  of  decreasing 
wavenumbers  k.  Asymptotic  analysis  shows  that  as  k  -+  0,  vortices  are 
confined  to  a  wall  layer  of  thickness  and  Vb  tends  to  a  constant 

when  moving  away  from  the  wall  layer. 
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Fig.  5  The  eigenfunctions  of  (3.32)  corresponding  to  a  set  of  increasing 
wavenumbers  k.  Asymptotic  analysis  shows  that  as  k  — *  oo,  vortices  are 
confined  to  an  0(fc-1/2)  thin  layer  centred  at  t;  =  2.3228  where  -T' /T2 
attains  its  maximum. 


53 


Local  Gortler  number 


Fig.6.  Neutral  curves  obtained  when  we  impose  the  following  initial  coni- 
tions  at  the  upstream  postion  Xo  =  60: 

(J  =  e-1^"  V'  =  0,  0  =  r/*7. 


Local  Gortler  number 


•  •  Asymptotic  result 

•  •  •  The  third  initial  condition 
— •  The  second  initial  condition 
—  The  first  initial  condition 
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Fig. 7.  Neutral  curves  obtained  when  we  impose  the  following  three  initial 
conditions  at  the  downstream  location  .Vo  =  20  and  use  £s  as  the  energy 
measure  to  monitor  the  energy  growth. 

(i) .  U  =  e-1/”"  •  e-'/r),  V  =  0,  9  =  r,U. 

(ii) .  U  =  T)2e-10/r,i  •  e-*»,  V  =  e"10^5  •  e"\  9  =  0, 

(iii) .  0  =  e~40^6  ■  e~n ,  V'  =  r?3e-40^6  •  e-r?s .  9  =  0. 
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Fig. 8.  Neutral  curves  obtained  when  we  impose  the  following  initial  con¬ 
ditions  at  three  different  upstream  positions  A'o  =  20,40,60  and  use  £5  as 
energy  measure: 

V  =  e~1/,7,<  • e~'n/r V  =  0.  8  =  tjU . 
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the  wall  mode 


Fig.  13  The  neutral  curves  for 
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Fig.  19  Variation  of  B  with  respect  to  the  wall  cooling  coefficient  n 
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FIg.20  A  comparison  of  neutral  curves  corresponding  to  two  models 

(El  as  energy  measure) 
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Fig.21  A  comparison  of  neutral  curves  corresponding  to  two  models 

(E2  as  energy  measure) 
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(E3  as  energy  measure) 
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local  Gortler  number 


Fig.  24  The  neutral  curve  for  Gortler  vortices  in  hypersonic  boundary  layer*.  The  dashed  part  of  the  curve 
indicates  that  at  the  given  wavenumber  the  curve  depends  at  zeroth  order  on  nonparalle  effects.  The  verical 
artows  indicate  the  size  of  the  next  order  correction  to  the  neutral  curve,  apart  from  the  regime  where  the 
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temperature  deeream.  mouotouically  to  lU  free  stream  value.  Gortler  eordcee  which  have  wavelength  comparable  with  the  bouudary  l»y«r  thictoe 
locU  wavenumber  of  order  U~'n)  Me  referred  to  ue  wull  mode..  We  .how  that  their  downetnam  evolution  w  goren.ee)  by  .  eet  of  pwbolic  pnrtral  differential 
equation.  end  that  they  hare  the  anal  featnre.  of  Gortler  aortic™  in  iacospreasble  bonadary  layer..  A.  the  local  wnrennmber  uicreaw,  the  nentral  Gortler 
namber  decree™,  end  the  centre  of  eortex  aetirity  more.  towMdj  the  temperature  adjustment  layer.  Gortler  rorticM  with  warennmber  of  order  one  or  larger 
met  oecemarily  be  trapped  in  the  temperature  adjustment  layer  and  it  i.  thb  mode  which  ia  the  mom  dangerou..  For  that  mode,  we  (tad  that  the  lending  order 
term  in  the  GSrtler  number  expanwon  i.  utdependent  of  the  wa renumber  and  u  due  to  the  currature  of  the  bamc  .tale.  Thi.  term  i.  aln  the  a^mptoUc  limit  of 
the  neutral  Gortler  number,  of  the  wall  mode.  To  determine  the  higher  order  correction  term,  ia  the  Gortler  number  expaamon.  we  hare  to  dutinguiih  between 
two  wall  currature  emu  When  the  wall  currature  ia  proportional  to  (Jur*'1  where  <  ia  the  rtreamwme  rnrinble,  the  Mach  number  if  can  be  wiled  oil  of  the 
problem  and  we  dow  that  in  the  0(1)  wareuumber  regime,  Gortler  rorticn  are  again  gorerued  by  a  act  of  parabolic  partial  differential  “l»*bon.  and  therefore 
the  higher  order  correction  term,  in  the  Gortler  namber  expaamon  are  sot  uniquely  determined  and  are  strongly  dependent  os  nonparallel  effect..  In  the  large 
warennmber  limit,  howerer.  nonparallel  effect,  become  of  aeeond  order;  Gortler  rortiem  erolre  downstream  in  a  qnawpurallel  manner  «d  the  Gortler  .amber 
expaamon  has  iU  first  three  term.  independent  of  nonparallel  effect..  In  the  more  general  case  when  the  wall  currature  is  not  proportional  to  ( Je)  the  effect  of 

the  currature  of  the  basic  state  permit,  in  the  downstream  derelopment  of  Gortler  rordoe.;  aon-paruilel  effect,  are  important  oeer  a  larger  range  of  ware. umber, 
ud  they  become  a  Mcond  order  only  when  the  w.renamber  i.  of  order  higher  than  0(Jf In  the  latter  cum  the  Gortler  number  expanmon  h™  the  ffrat  two 
term,  independent  of  nonparallel  effect.;  the  first  term  being  doe  to  lie  currature  of  the  bade  .Ute  and  the  second  term  due  to  rueou.  effect*.  The  mcond  term 
become,  comparable  with  the  fir.,  term  when  the  w.renamber  reache.  the  order  V*'V  in  which  ram  another  correction  term  can  abo  be  toaad  independently  of 
nonparallel  effect.  Next  we  mrmtig.te  real  gu  effects  by  winning  th.t  the  fiuid  i.  as  ideal  dimod.Ung  g“  We  find  that  both  gu  dimociation  and  wMl  cooling 
are  dratabiluing  for  the  mode  trapped  in  the  temperatnre  adjustment  l.yer,  bat  tot  the  wall  mode  trapped  Bear  the  wall  the  effect  of  ga.  diMocrauoi  can  be 
either  dotsbilimsg  or  stohilimsg. 
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